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The  EPAEC  gas  washout  model  has  been  applied  to  the  organic  tracers  ethyl- 

i 

acetoacetate  and  diethylamine.  This  model  utilises  source  strength,  vapor 
solubility  and  basic  meteorological  parameters  to  estimate  the  tracer  cori- 
centration  in  downwind  rain  samples.  Contrary  to  the  study  of  SO2  washout, 
the  observed  concentrations  are  in  rather  poor  agreement  with  predictions 
derived  from  the  model.  A  detailed  error  analysis  has  identified  poorly 
defined  tracer  solubility  as  the  major  source  of  this  discrepancy,  v Solu¬ 
bilities  expressed  in  terms  of  the  Henry's  Law  relationship. 


H' 


vapor  pressure 


,  have  been  derived  from  the  available  literature 


molar  concentration 
and  are  judged  to  be  inappropriate  for  the  low  vapor  concentrations 
encountered.^  Field  results  have  also  been  used  to  calculate  the  fractional 
washout  of  the  vapor  per  kilometer  from  the  source.  These  results  together 
with  existing  diffusion  models  may  be  used  to  estimate  the  vapor  concentra¬ 
tions  of  these  materials  along  a  plume  trajectory. 
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Th*  pr 'miry  objective 


Ls  research  program  is  to  evaluate  an- 


procipi  s,1t  ion  scavenging  as  a  mechanism  for  removal  of  tracp  organ  iC 

^  il  ... 


gases  ar.d  vapors  from  the  atmosphere.  Significant  levels  of  organic 
vapors  in  various  stages  of  oxidation  exist  in  both  polluted  and  natural 
atmospheres,  and  the  existence  of  measurable  amounts  of  organic  materials 
in  natural  rainwater  has  been  recognized  for  some  time.  Junge  (1963), 
for  instance,  indicates  that  about  1  milligram  of  organic  per  liter  of 
rain  may  be  expected  under  average  rural  pollution-level  and  rainfall 
conditions.  This  would  mean  that  for  over  a  land  surface  of  one  square 
kilometer,  each  millimeter  of  rain  scavenges,  under  average  rural  con¬ 
ditions,  about  one  kilogram  of  organic  material. 


Presumably  the  degree  of  organic  scavenging  varies  with  the  concentration 
and  type  of  material  in  question,  although  practicallv  no  research  has 
been  dedicated  to  this  phenomenon  to  date  insofar  as  gases  and  vapors  are 
concerned.  Also  very  little  attempt  has  been  made  to  isolate  specific 
organics  within  natural  rainfall  samples,  although  the  limited  work  tnat 


has  been  done  in  this  area  interestingly  indicates  that  along  with  the 


simpler  organics  some  surprisingly  complex  forms  exist,  including  poly¬ 
saccharides,  organochlorine  compounds,  and  vitamin  B ^ 2  (cf*  Semenov,  - 
al  (1967)  (1968),  Tarrant  and  Tatton  (1968),  and  Parker  (1968)). 


Significant  experimental  research  lias  been  carried  out  regarding  the 
scavenging  of  organic  particulates.  These  aerosols,  however,  which  have 
been  confined  primarily  to  dye  materials  because  of  their  utilitv  as 
tracers,  are  not  expected  to  show  significantly  different  behavior  from  inor¬ 
ganic  aerosols  having  similar  physical  characteristics.  Because  of  this, 
the  ongoing  and  past  research  involving  aerosol  scavenging  can  be  con 
sidered  to  apply  equally  well  to  both  organics  and  inorganics,  and  at 


G 


present,  a  comparatively  broad  base  of  scientific  information  exists  in 


this  area. 


This  same  statement  cannot  be  made  with  regard  to  organic  gases  and 
vapors,  however.  Past  studies  of  gas  scavenging  have  been  limited  to 
inorganic  materials,  and  only  recently  has  their  level  of  effort  been 
advanced  to  a  state  comparable  with  that  of  aerosol  research.  During  the 
past  ten  years,  however,  significant  and  encouraging  advances  in  the 


field  of  gas  scavenging  have  been  made,  including  the  field  studies  of  I2 
and  Br2  washout  by  Engelmann  (1968),  the  laboratory  investigations  of 
SO2  and  NC>2  washout  by  Georgi  and  Bielke  (1970),  the  theoretical  analysis 
of  tritium  washout  by  Hales  (1972a)  and  the  field  and  theoretical  anal¬ 
ysis  of  SO2  washout  by  Hales  and  his  coworkers  (1972b)  (1972)  (1973). 


These  latter  workers  pointed  out  and  experimentally  verified  the  fact 
that  for  gas  scavenging,  both  absorption  of  the  gas  into  the  raindrops 
and  its  desorption  from  the  raindrops  are  important  in  the  overall  scav¬ 
enging  process.  This  "reversible"  aspect  depends  primarily  upon  the 
solubility  of  the  gas  and  contrasts  with  the  behavior  of  aerosols,  which 
can  reasonably  be  assumed  to  attach  permanently  to  raindrops  in  an 
"irreversible"  manner. 


Demonstration  of  the  reversible  interaction  between  trace  gases  and 
natural  rain  led  to  the  development  of  a  descriptive  model  for  gas  scav¬ 
enging.  Although  generalize* ' e  to  other  source  configurations,  this  model 
envisioned  at  the  outset  a  plume  emitted  from  a  single  point  source  as 
shown  in  Figure  1.  Rain  falling  through  the  plume  absorbed  gas  as  it 
fell  into  the  concentrated  regions,  and  under  some  conditions  desorbed 
material  as  it  fell  through  the  less  concentrated  areas  below.  The 
theoretical  treatment  was  based  on  conventional  mass-transfer  theory; 


for  the  set  of  simplifying  conditions  including 

1.  A  Gaussian  plume  with  invariant  background  concentration 


2.  Vertical  rainfall 


3.  Linear  mass  transfer  and  solubility  behavior 


>  vw*  V*» 


o 


4.  First-order  gas-phase  reaction  kinetics, 
this  model  could  be  expressed  explicitly  in  terms  of  the  equations 


3QK  r  /  2  o2i2 

'  ■  2 exp  r  ^ + -5“ 

2V2ti  o  uv  a  V  y 

y  z  \ 


exp  (ch)  If  1-erf 


+  exp  (-  Ch)  1-erf 


+  WH' 


3k  H' 
_  X_ 

av 

z 


~v~h  \ 

a  vr/ 


-°2^ 


r  *  exp  (-kx/u) . (3) 


where  c^(a)  Is  the  concentration  of  pollutant  in  radius-a  drops  at  the 
ground-level  receptor  point.  Other  terms  in  the  equation  are  defined  in 
the  Table  of  Nomenclature  and  are  summarized  here  as  follows: 

h  *  release  height 

H'  =  solubility  constant 

k  *  reaction  rate  constant 

Ky  *  overall  mass-transfer  coefficient 

Q  *  emission  source  strength 

r  ■  rate  of  chemical  reaction  of  pollutant 

vz  =  terminal  fall  velocity  of  raindrop 

y^k  ■  background  pollutant  mixing  ratio 

o„,  o_  *  plume  dispersion  parameters 
y  z 
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A  more  accurate  version  of  this  model  was  subsequently  developed,  which 
was  not  constrained  by  assumptions  2  and  3  above.  This  version,  which 
has  been  named  the  "EPAEC"  scavenging  model,  cannot  be  presented  in  an 
explicit  form  similar  to  that  of  Equation  1.  It  can  be  described  super¬ 
ficially,  however,  in  terms  of  the  computer  subroutine  hierarchy  shown  in 
Figure  2.  This  algorithm  computes  ground  level  concentrations  numerically 
comparable  to  those  calculated  analytically  by  Equation  1.  It  also  pro¬ 
ceeds  further  to  calculate  weighted  averages  of  concentrations  in  accord¬ 
ance  with  prescribed  raindrop  size  statistics,  and  computes  cumulative 
crosswind  washout  values.  A  feature  of  the  EPAEC  model  is  its  modular 
construction,  which  enables  rapid  modification  and  maximum  utility.  Thus, 
if  calculations  for  a  different  type  of  plume  (a  puff,  say)  were  desired, 
one  would  simply  replace  the  plume  model  module  and  proceed  to  execute 
the  program.  Similarly  a  variety  of  gases  can  be  accommodated  by  utiliz¬ 
ation  of  the  appropriate  solubility  model  in  the  overall  program.  A 
detailed  description  of  the  EPAEC  model  has  been  presented  in  an  earlier 
report  (Dana  (1973))  and  is  included  with  this  document  as  an  appendix. 

SCIENTIFIC  APPROACH 

The  EPAEC  model  is  considered  to  be  a  rather  realistic  description  of 
precipitation  scavenging  of  gases  from  plumes.  It  has  been  shown  to  com¬ 
pare  favorably  with  experimental  measu  ;.ients  of  SO2  washout  under  a 
variety  of  conditions,  indicating  its  potential  utility  as  a  means  of 
investigating  (and  subsequently  modeling)  the  washout  of  organic  mater¬ 
ials.  The  scientific  approach  of  this  program,  therefore,  is  to  perform 
field  experiments  of  precipitation  scavenging  of  organic  gases  under 
well-defined  conditions  and  compare  these  results  with  corresponding 
predictions  of  the  EPAEC  model.  Deviations  between  experiment  and  theory 
then  an  be  analyzed  in  view  of  basic  model  assumptions  and  subsequent 
improvements  can  be  attempted  on  the  basis  of  this  analysis.  The 
following  chapter  presents  a  detailed  description  of  these  comparisons. 
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Figure  2.  SUBROUTINE  HIERARCHY  IN  EPAEC  MODEL. 
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CHAPTER  II. 

FIELD  EXPERIMENTS 


Tne  tests  described  in  this  report  were  conducted  at  the  ''uillavute  Air 
field  on  the  Olympic  Peninsula  of  western  Washington.  This  site  offered 
the  advantages  of  ample  rainfall,  a  U.  S.  Weather  Bureau  station  and 
isolation  from  pollution  sources  which  might  contribute  significant  back¬ 
ground  interference.  While  much  of  the  terrain  in  this  region  is  heavily 
wooded,  the  airfield  provided  a  sizable  area  free  of  structures  which 
would  yield  irregularities  in  the  air  motion  over  the  sampling  network. 
Although  hangars  located  approximately  a  quarter  mile  upwind  of  the  tower 
used  to  release  the  tracers  mav  have  been  the  source  of  turbulence  observ¬ 
ed  during  previous  tests  at  this  site,  meteorological  data  recorded 
during  this  series  provided  no  evidence  of  such  turbulence.  The  stand-bv 
status  of  the  airfield  allowed  tests  to  be  conducted  without  interference 
during  daylight  hours.  However,  its  designation  as  an  emergency  landing 
strip  required  the  tower  to  be  lowered  during  the  evening  hours  in  order 
to  eliminate  it  as  a  hazard  to  aircraft  requiring  use  of  the  airfield 
during  this  period  of  reduced  visibility.  This  stipulation  limited  the 
number  of  tests  which  could  be  conducted  and  extended  tie  time  required 
to  put  the  release  equipment  into  operation. 

Figure  3  illustrates  the  distribution  of  the  rain  samplers  ith  respect 
to  the  tower.  Arc  A,  at  300  feet,  consisted  of  17  sampling  sites  located 
at  5°  intervals  within  the  80°  arc  defined  bv  lines  drawn  NNE(20°)  and 
WNW(300°)  from  the  tower.  Arc  B  consisted  of  29  sampling  locations  en¬ 
compassing  the  same  angle  as  arc  A  and  following  the  geometry  of  the 
airfield  runways.  Although  the  prevailing  winds  during  each  test  carried 
the  tracer  through  the  region  described  bv  the  arcs,  additional  samnlers 
were  set  out  during  experiments  2,  3,  4  and  5  to  insure  nlume  enclosur**. 
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Rain  samples  were  collected  in  polyethylene  bottles  fitted  with  7  3/8- 
inch  polyethylene  funnels  as  shown  in  Figure  4.  In  order  to  inhibit 
desorption  of  the  volatile  tracer,  the  rain  was  frozen  as  it  was  collected. 
This  was  accomplished  by  thermally  insulating  the  sample  bottles  in  polv- 
styrene  containers  and  surrounding  them  with  dry  ice. 

In  addition  to  the  collection  of  rain  samples,  the  integrated  vapor 
phase  concentration  of  tracer  at  ground  level  was  determined  through  the 
use  of  a  limited  number  of  vapor  samplers,  Figure  4.  These  units  con¬ 
sisted  of  a  battery-operated  pump  which  drew  air  at  1  1pm  through  an 
impinger  containing  10  ml  of  distilled  water  (ethylacetoacetate  release) 
or  10  ml  of  5  x  10“^  N  nitic  acid  (diethylamine  release).  At  the  concl  i- 
sion  of  an  experiment,  the  impinger  solutions  were  transferred  to  polv- 
ethylene  bottles  for  storage.  The  impingers  were  then  rinsed  with 
distilled  water  in  preparation  for  the  next  test.  Addition  of  fresh 
distilled  water  or  dilute  acid  was  delayed  until  the  period  immediately 
preceding  the  next  experiment.  Analysis  of  the  impinger  solution  gave 
the  total  mass  of  tracer  in  the  air  sampled  and  provided  an  indication  of 
the  contribution  of  dry  deposition  to  the  tracer  concentration  of  the  rain¬ 
water.  Since  only  a  limited  number  of  the  vapor  samplers  were  available, 
their  distribution  among  the  sampling  sites  was  dictated  by  the  observed 
wind  conditions  prior  to  release  of  the  tracer. 

Of  major  concern  in  determining  the  success  of  the  washout  experiment  was 
the  tracer  dispersion  system.  Since  one  of  the  tracers  was  not  especiallv 
volatile  (boiling  point  180°),  it  was  essential  that  complete  evaporation 
be  assured  while  maintaining  a  rate  of  release  which  was  sufficient  to 
allow  detection  of  the  tracer  in  the  downwind  precipitation  samples. 

The  tracer  dispersion  equipment  designed  and  built  by  Batcelle-North- 
west  to  meet  the  above  requirements  is  shown  in  Figures  5  and  6.  It 
consists  of  a  5 -gallon  storage  tank,  a  com]rr's>ed  air  cylinder  and  vanor 
gun  containing  two  sonic  nozzles  mounted  ins-  .  le  an  8-inch  aluminum  tube. 
Release  was  made  by  transferring  the  tracer  under  low  pressure  to  the 
sonic  nozzles  and  injection  into  the  atmosphere  bv  high  pressure  air 
from  an  air  compressor.  A  heat  gun  was  positioned  directly  behind  the 
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Figure  4.  RAINWATER  AND  VAPOR  SAMPLERS 


Figure  5. 


VAPOR  DISPERSION  ASSEMBLY:  1.  TRACER  STORAGE 
TANK;  2.  NITROGEN  FOR  LOW  PRESSURE  TRANSFER  OF 
TRACER; 3.  AIR  COMPRESSOR;  4.  VAPOR  GUN. 
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sonic  noczles  so  thst  the  sir  stresm  receiving  the  tracer  was  slightly 
shove  ambient  temperature.  Ground  level  testing  of  this  equipment  is 
shown  In  Figure  7. 

When  favorable  test  conditions  were  observed,  the  tower  was  raised,  the 
tracer  storage  tank  filled  and  the  vapor  dispersion  equipment  secured  at 
the  top  of  the  tower  (100  ft) .  The  configuration  of  the  meteorological 
sensors  and  the  release  equipment  is  shown  in  Figure  8.  While  the  tower 
was  being  placed  in  position,  the  rain  and  vapor  samplers  were  distributed 
as  illustrated  previously.  In  order  to  confine  sampling  as  nearly  as 
possible  to  the  period  of  release,  the  implngers  were  turned  off  and  the 
rain  samplers  covered  until  tower  preparation  was  complete.  The  sampling 
equipment  was  then  placed  in  operation  before  generation  of  the  organic 
plume.  Although  some  rain  was  necessarily  collected  before  dispersion 
of  the  tracer,  the  effect  of  the  resulting  dilution  was  at  least  partially 
compensated  for  by  the  procedure  followed  in  retrieving  the  samples. 
Samples  were  collected  in  a  sequence  which  was  the  same  as  that  followed 
in  putting  them  into  service. 

In  addition  to  the  downwind  sampling  locations  identified  in  Figure  3, 
at  least  one  background  rain  and  vapor  sample  was  obtained  at  a  site 
>100  feet  upwind  of  the  source.  These  samples  were  handled  in  the  same 
manner  as  those  located  downwind  of  the  tower.  The  results  of  the 
chemical  analysis  of  all  field  samples  are  recorded  in  Appendix  D.  It 
will  be  noted  that  most  background  samples  appear  to  indicate  the 
presence  of  some  tracer.  The  origin  of  this  apparent  tracer  load  is 
discussed  In  Chapter  5.  It  is  evident,  however,  that  the  magnitude 
of  this  background  interference  is  too  low  to  obscure  the  washout 
behavior  of  the  tracers. 
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CHEMICAL  ANALYSIS 

ANALYTICAL  PROCEDURES 

Analysis  of  field  samples  was  accomplished  by  colorimetric  methods  which 
are  documented  in  the  chemical  literature.  The  initial  phase  of  this 
study  was  directed  to  the  investigation  of  these  procedures  in  order  to 
determine  their  respective  detection  limits  and  identify  problems  which 
might  arise  from  the  methods  to  be  used  for  the  collection  and  storage  of 
the  field  samples.  This  preliminary  work  as  well  as  the  subsequent  analy¬ 
sis  of  field  samples  was  done  with  a  Bausch  and  Lomb  Spectronic  20 
spectrometer  using  a  1-inch  test  tube  as  the  sample  cell.  A  discussion 
of  this  preliminary  work  will  provide  an  indication  of  the  confidence  one 
can  place  in  the  results  tabulated  in  Chapter  5  and  follows  a  description 
of  the  ethylacetoacetate  and  diethylamlne  analysis  procedures. 

Ethylacetoacetate  Analysis — A  25.0  ml  aliquot  of  the  ethylacetoacetate 
solution  was  pipetted  into  the  spectrometer  cell  followed  by  2.0  ml  of 
chilled  dlazonium  salt  solution  prepared  by  mixing  equal  volumes  of  0.5Z 
potassium  nitrite  and  0.007A  M  p-aminoacetophenone  (0.1  M  in  HC1).  Addi¬ 
tion  of  1.0  ml  of  reagent  grade  piperidine  to  the  well-mixed  solution 
produced  a  transient  dlazonium  addition  product  having  an  absorption 
maximum  at  A35  nm.  Observation  of  the  change  in  transmittance  with  time 
confirms  the  report  published  in  the  Operation  Procedure  Guide  (1965)  that 
maximum  color  development  occurs  3  minutes  after  piperidine  introduction, 
and  does  not  undergo  noticeable  decay  for  a  period  of  more  than  1  minute. 
Transmittance  readings  were  therefore  taken  3  to  A  minutes  after  piperi¬ 
dine  addition. 

The  analysis  of  standard  ethylacetoacetate  solutions  shows  the  Beer's  Law 
plot  to  be  linear  to  2A  yg  (0.96  ppm)  with  a  detection  limit  of  1.0  yg 
(0.0*  ppm).  Analysis  of  those  solutions  which  show  concentrations  greater 
than  0.9  ppm  were  repeated  using  a  smaller  sample  aliquot,  diluting  to 
25.0  ml  and  analyzing  as  before. 

*  An  estimate  of  the  analytical  error  was  obtained  from  analysis  of  stand¬ 
ard  solutions  included  with  each  batch  of  field  samples.  These  solutions, 
spanning  the  linear  Beer's  Law  range,  show  a  standard  deviation  of  less 
than  25  per  cent. 
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Dlethylamine  Analysis — The  colorimetric  procedure  reported  by  Dahrgren 
(1964)  for  the  analysis  of  aqueous  amine  solutions  was  adopted  for  this 
study.  A  25.0  ml  aliquot  of  the  sample  solution  was  pipetted  into  the 
spectrometer  cell  and  the  pH  adjusted  by  the  addition  of  1.0  ml  of  phos¬ 
phate  buffer  (pH  6.6).  One  milliliter  of  10“  ^  M  sodium  hypochlorite  was 
added  allowing  10  minutes  for  complete  conversion  of  the  .-.mine  before 
destroying  the  excess  hypochlorite  with  1.0  ml  of  0.5%  sodium  nitrite. 
Reaction  of  the  chloramine  product  with  1.0  ml  of  starch-potassium  iodide 
solution  yielded  the  intense  blue  color  characteristic  of  the  starch- 
triiodide  complex.  Fifteen  minutes  was  allowed  for  maximum  color  develop¬ 
ment  before  measuring  the  transmittance  at  540  nm. 

Examination  of  a  series  of  standard  diethylamine  solutions  shows  a  linear 
Beer's  Law  plot  to  at  least  10.0  ug  (0.5  ppm)  with  a  detection  limit  of 
approximately  2  ug  (0.08  ppm).* 

Since  all  field  samples  were  frozen  during  collection  and  stored  at  sub¬ 
freezing  temperatures,  these  solutions  were  allowed  to  warm  slowly  to 
room  temperature  before  proceeding  with  the  analysis  as  outlined  above. 

STABILITY  OF  AQUEOUS  TRACER  SOLUTIONS 

Among  the  questions  which  required  examination  before  proceeding  with 
the  field  tests  was  the  problem  of  tracer  stability.  As  with  all  experi¬ 
ments  dependent  on  accurate  chemical  analysis  of  environmental  samples, 
it  was  exceedingly  important  to  consider  those  mechanisms  which  might 
diminish  tracer  concentration  during  storage.  Since  the  rain  samples 
were  frozen  as  they  were  collected,  desorption  of  the  volatile  organic 
during  the  field  tests  dees  not  appear  to  be  a  viable  path  for  tracer 
loss.  However,  chemical  reaction,  such  as  hydrolysis  of  ethylacetoacetate, 

and  adsorption  to  container  walls  during  storage  were  recognized  as 

_ 

As  with  the  ethylacetoacetate  samples,  solutions  containing  known  amounts 
of  diethylamine  were  analyzed  with  the  field  samples.  This  served  as  a 
means  of  insuring  against  possible  positive  interference  by  ammonia  in  the 
reagent  solutions  and  allowed  an  estimate  of  precision.  These  results,  as 
well  as  multiple  analysis  of  some  field  samples,  suggest  a  maximum  exper¬ 
imental  error  of  approximately  30  per  cent. 


potentially  si leant .  Some  insight  into  the  role  of  these  processes 
was  gained  oy  preparing  dilute  solutions  of  the  tracers,  and  storing  them 
at  room  temperature  in  polyethylene  bottles  of  the  type  used  in  the  field 
experiments.  Analysis  performed  on  aliquots  of  these  solutions  were  in 
excellent  agreement  with  those  conducted  immediately  after  their  prepara¬ 
tion.  Some  of  the  data  acquired  during  this  preliminary  study  are  summ¬ 
arized  in  Table  1.  It  is  evident  from  this  work  that  dilute  tracer  solu¬ 
tions  are  sufficiently  stable  to  allow  them  to  be  stored  for  several  weeks. 
Storage  at  sub-freezing  temperatures  and  analysis  at  the  Bat  telle  Hanford 
laboratory  was  therefore  adopted  as  the  procedure  to  be  followed  in 
handling  the  field  samples. 


TABLE  1. 

STABILITY  OF  AQUEOUS  TRACER  SOLUTIONS 
Initial 

Tracer  Concentration  Concentration  (Days  of  Storage) 


Ethylacetoacetate 
E t hy lace toace tat e 
Diethylamine 


0.81  ppm 

0.81 

0.40 

0.43 

0.25 

0.26 

(4)  0.83  (24) 
(4)  0.42  (24) 
(29) 


The  vapor  sampling  systems  used  to  provide  a  measure  of  the  integrated 
plume  concentration  at  ground  level,  while  not  essential  to  the  field 
program,  did  provide  input  for  accessing  the  contribution  of  dry  deposi¬ 
tion  to  the  tracer  load  of  the  collected  rain  samples.  Since  the  utility 
of  these  systems  depended  on  the  efficient  removal  of  the  tracer  from  air 
drawn  through  a  midget  impinger,  preliminary  work  included  an  investiga¬ 
tion  of  these  systems.  These  tests  consisted  of  placing  a  dilute  solution 
of  the  tracer  in  the  impinger  and  drawing  clean  air  through  it  for  a 
period  of  approximately  two  hours.  The  impinger  solutions  were  analyzed 
to  determine  the  evaporative  loss  of  tracer  during  operation  of  the  sampler. 
As  with  the  field  experiments,  distilled  water  was  used  for  the  ethylaceto¬ 
acetate  test  and  dilute  nitric  acid  (5  x  10“^  N)  for  the  diethylamine  test. 
The  results  recorded  in  Table  2  illustrate  that  dilute  acid  forms  an 
efficient  trap  for  diethylamine  capture  while  distilled  water  is  somewhat 
less  effective  in  scrubbing  ethylacetoacetate  from  air  drawn  through  the 


0 


imp Inge r.  Samples  obtained  during  Experiments  1  and  4  have  been  analyzed 
and  the  calculated  mass  of  tracer  In  each  sample  recorded  In  Tables  14 
and  18. 
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TABLE  2. 

THE  RESULTS  OF  TESTS  DESIGNED  TO  INDICATE  THE  EFFICIENCY 
OF  THE  VAPOR, SAMPLING  SYSTEM 


o 


0 


Implnger  yg  Tracer  Period  of 


Tracer  Solution  Initial 

Final 

Operation 

Per  Cent  Loss 

distilled  water 

Ethylacetoacetate  20 . 2 

12.3 

2  hr.  15  min. 

39 

5  x  10“5  N  HN03 

Dlethylamlne  5.0 

5.1 

2  hr. 

0 

O 


O 


O 


O 
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CHAPTER  IV. 

DATA  PROCESSING 

As  indicated  In  Chapter  II,  a  number  of  support  measurements  were  per¬ 
formed  in  addition  to  the  assessment  of  washout  concentrations.  These 
included  the  anemometer  measurements,  rain-gauge  records  and  raindrop-size 
spectrum  determinations.  This  section  describes  the  procedures  used  for 
processing  these  data. 

The  Gill  anemometer  data  were  processed  in  several  stages.  The  first  of 
these  steps  involved  reading  the  field-generated  METRODATA  tapes,  conver¬ 
ting  the  raw  data,  and  writing  on  industry-compatible  7-track  tapes  using 
a  Bat  tell.;  Atmospheric  Sciences  Department-owned  NOVA  computer.  This 
system  employed  a  data  culling  routine  which  detected  and  flagged  faulty 
data  entries  on  the  output  tapes,  so  that  tape  quality  could  be  monitored 
visually  by  observing  printed  listings. 

The  7-track  tapes  were  listed  using  a  high  speed  printer  in  conjunction 
with  an  AEC-owned  UNIVAC  1108  computer.  Thesp  listings,  which  were  in 
the  form  of  tables  giving  clock  readings  and  uncalibrated  values  of  the 
U,  V,  and  W  wind  components  at  0.42  second  intervals,  were  scrutinized 
for  faulty  data  as  detected  by  the  NOVA  culling  routine.  These  checks 
showed  that  the  wind  data  had  been  recorded  with  high  quality;  the  few 
culled  data  points  that  did  appear  in  the  records  were  replaced  with  the 
values  immediate  preceding,  resulting  in  essentially  zero  distortion 
of  the  processed  results. 

Following  the  quality-control  checks,  the  7-track  data  tapes  were  resub¬ 
mitted  to  the  UNIVAC  computer  where  calibrated  mean  wind  velocities, 
directions,  and  standard  deviations  of  vertical  and  eziruthal  wind  direc¬ 
tions  were  calculated.  Computed  standard  deviations  were  based  on  sampling 
and  averaging  times  of,  respectively,  the  experiment  times  and  one-fourth 
the  source- ieceptor  transit  times.  Formulae  employed  for  this  purpose 
are  listed  as  follows: 
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The  nomenclature  in  these  equations  is  consistent  with  that  used  through- 
ou:  the  meteorological  literature  and  is  defined  in  the  Table  of  Nomen¬ 
clature.  Results  of  these  calculations  for  the  anemometer  data  taken 
from  the  field  experiments  of  this  study  are  given  in  Table  3. 


TABLE  3. 

METEOROLOGICAL  PARAMETERS  CALCULATED  FROM  GILL  ANEMOMETER  DATA 


Run 

Sampling 

Time 

Sec. 

Average 
Wind 
Speed 
cm/s ec 

Average 

Heading, 

(From  True  North) 
Radians 

°e 

Radians 

Radians 

1A* 

300 

A16.6 

2.68 

0.159 

O.OA35 

IB 

1210 

A67.2 

2.79 

0.22A 

0.0727 

2 

1A90 

328. A 

2.35 

0.233 

0.067A 

3 

2390 

376.3 

2.93 

0.137 

0.0757 

A 

1820 

711.2 

3.00 

0.09A2 

0.06A8 

5 

1620 

711.7 

3.00 

0.105 

0.0676 

*Run  1  divided  into  two  components  owing  to  temporary  lull  in  rain  and 
consequential  shutdown  during  mid-experiment. 

Operation  of  the  raindrop  spectrometer  provided  ozalid  images  as  shown  hv 
the  example  in  Figure  9.  The  sizes  of  these  raindrop  images  were  related 
to  actual  drop  sizes  by  a  previous  calibration:  thus  the  raindrop  spectra 
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Figure  9.  REPRESENTATIVE  OZALID  IMAGES  FROM  (a)  EXPERIMENT  3  AND 
(b)  EXPERIMENT  5. 
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could  be  obtained  simply  by  sizing  the  Images  and  applying  the  calibra¬ 
tion.  Image  sizing  was  performed  using  a  Zeiss  Spectrometer.  Resulting 
raindrop  size  spectra  determined  for  each  of  tin  five  tests  are  presented 
in  Table  4. 


0 


Table  4.  MEASURED  RAINDROP  SIZE  SPECTRA  AND  RAIN  RATES. 


Size  Range  (Diameter) 


mm 

Run  1 

Run  2 

Run  3 

Run  4 

Run  5 

0 

0- 

224 

0.122 

0.305 

0.134 

0.170 

0.237 

.224- 

.349 

0.268 

0.389 

0.230 

0.285 

0.180 

.349- 

.417 

0.069 

0.058 

0.111 

0.090 

0.094 

.417- 

.497 

0.147 

0.033 

0.065 

0.050 

0.075 

0 

.497- 

.594 

0.182 

0.028 

0.054 

0.035 

0.043 

.954- 

.710 

0.081 

0.033 

0.122 

0.065 

0.113 

o 

1 

.847 

0.039 

0.051 

0.094 

0.110 

0.075 

1 

*3" 

00 

1.011 

0.036 

0.045 

0.088 

0.085 

0.056 

• 

1.011- 

1.321 

0.042 

0.051 

0.074 

0.070 

0.082 

1.321- 

2.117 

0.009 

0.007 

0.028 

0.040 

0.045 

Mass  Mean 
Drop  Size 


• 

mm 

0.602 

0.590 

0.726 

0.740 

0.772 

Rain  Rate 
cm  /hr 

0.20 

0.59 

0.25 

0.25 

1.1 

In  addition  to  showing  the  raindrop  spectra,  Table  4  presents  the  mean 
rainfall  rates  recorded  onsite  by  the  fast-response  rain  gauge.  The  data 
obtained  from  this  unit  were  intended  primarily  for  the  purpose  of  analvz- 
ing  the  effects  caused  by  the  finite  time  intervals  required  for  deploy¬ 
ment  of  the  sampling  arcs:  actual  charts  from  the  rain  gauge  for  each  test 
are  reproduced  in  Appendix  C. 


Computations  with  the  F.PAEC  model  required,  in  addition  to  the  data 
provided  in  Tables  3  and  4,  information  pertaining  to  the  molecular 
transport  and  solubility  properties.  Values  of  the  specific  properties 
needed  for  this  purpose  are  given  in  Table  5. 


I 
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Table  5.  SOLUBILITY  AND  MOl.PCtfl.AJ*  TRANSPORT  PROPERTIES. 


Ethylacetoacetate 

Diethylamine 

Diffusivity  in  Air 
cm^/sec 

.074 

.096 

Solubility  in  Water* 
at  10°  C 

moles /cm^ 

4.0 

.091 

Kinematic  Viscosity 
of  Air 
cm^/sec 

.14 

1  c 

1  c 

*  Defined  by:  Solubility  M  — »  where  c  is  the  liquid 

phase  concentration  in  mcles/cm^  and  y  is  the  gas-phase 
mole  fraction. 


Diffusion  coefficients  of  the  test  materials  were  calculated  using 
collision  integral  theory;  associated  errors  are  expected  to  be  of  the 
order  of  five  percent.  Solubilities  appearing  in  Table  5  are  known  with 
much  less  certainty,  perhaps  being  in  error  by  factors  of  five  or  more. 

The  reason  for  this  uncertainty  is  the  unanticipated  lack  of  solubility 
data  reported  for  these  materials  at  the  low  concentrations  of  present 
interest.  These  values  are  "best  estimates"  obtained  from  high-concentra¬ 
tion  data  by  assuming  the  validity  of  Henry's  Law  and  applying  the  Clausius 
Claperyon  Equation.  A  detailed  description  of  the  calculations  leading  to 
the  values  in  Table  5  is  given  in  Appendix  A. 

A  description  of  the  actual  entry  of  data  into  the  EPAEC  model  can  be 
facilitated  by  considering  the  input  segment  of  the  program  reproduced 
below: 


9 
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100 

READ 

(5,320) 

N, J1,J2,J3,J4,  JOPT, JP, JEND 

IF 

(Jl.EQ.l) 

GO  TO  110 

READ 

(5,330) 

(D(I),I-1,  N) 

READ 

(5,330) 

(F(I),I»1,  N) 

110 

IF 

(J2.EQ.1) 

GO  TO  120 

READ 

(5,340) 

DAX , DAY ,  HEX , XN U , P , T , XNT , RK 

120 

IF 

(J3.EQ.1) 

GO  TO  130 

READ 

(5,350) 

SIGTHE,  SIGPHI, U.H.Q, VERT, BKG 

130 

IF 

(J4.EQ.1) 

GO  TO  140 

READ 

(5,360) 

XBUK, YBUK, ZBUK,  DELTAY ,  DELTAZ 

o 


o 


o 


0 


o 


The  first  READ  statement  reads  the  number  of  steps  in  the  discrete  rain¬ 
drop  size  spectrum  (N)  in  addition  to  a  series  of  control  variables. 

If  these  are  set  as  follows: 

J1-J4  =  0 

JOPT  *  1  for  gas-phase  limited  conditions.  JOPT  =  0  for 
stagnant-drop  conditions. 

JP  -  1 
JEND  «  0. 

The  variables  used  in  this  code  are  described  in  detail  in  Appendix  E 
and  will  only  be  summarized  here,  as  follows: 

D,  F  -  Raindrop  Diameter  and  Frequency 

DAX*,  DAY  -  Liquid  and  Gas  Phase  Diffusivities 

HEX*  -  Solubility  Parameter 

P,  T,  XNT  -  Pressure,  Temperature,  Rain  Rate 

RK*  -  Reaction  Rate  Constant 

SIGPHI ,  SIGTHE,  U,  H,  Q  -  Plume  Parameters 

VERT*  -  Plume  Loft  Velocity 

BKG*  -  Background  Concentration 

XBUK,  YBUK,  ZBUK,  DELTAY,  DELTAZ  -  Receptor  Location  and  Grid 

Spacing  Parameters 

A  typical  listing  of  input  data  is  given  below: 


5  0 

0  0 

0  0 

1  1 

.01 

.02 

.03 

.04 

.05 

.05 

.10 

.25 

.45 

.15 

. 74E-1 

1.4E-1 

1.0EO 

.062 

.056 

5.13E2 

3.05E3 

6.6E-2 

3.0E4 


Variable  not  significant  to  present  useage  of  program. 
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Solubility  data  were  entered  i..  subroutine  form  rather  than  through  READ 
®  statements.  The  EPAEC  routine  was  executed  several  times  for  each  run  to 

test  the  effects  of  varying  solubilities.  These  results  are  presented 
for  each  field  experiment  in  the  following  chapter. 

o 


o 


o 


o 


o 


o 


0 
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CHAPTER  V. 


RESULTS 

In  the  preceding  chapters  we  have  attempted  to  place  this  work  with 
respect  to  previous  investigations  of  washout  and  the  contribution  of 
this  study  to  the  development  of  a  comprehensive  washout  model.  We  have 
also  sought  to  identify  and  explore  questions  having  a  direct  bearing  on 
the  results  of  this  study  and  hence  the  interpretations  and  conclusions 
to  be  drawn  from  it.  In  this  chapter  we  will  discuss  the  five  field  ex¬ 
periments  which  focus  on  the  washout  behavior  of  ethylacetoacetate  and 
diethylamine  vapors.  A  description  of  the  meso  and  micro  scale  neteorologv 
will  be  followed  by  descriptive  observations  taken  from  the  field  notes 
and  a  graphical  presentation  of  the  analytical  data.  The  observed  tracer 
concentrations  in  rain  water  will  then  be  compared  with  those  predicted  by 
the  washout  model.  Since  the  tracer  solubilities  calculated  in  Appendix  A 
represent  only  estimates  based  on  the  best  available  literature  data, 
concentrates  are  also  preo^cted  assuming  a  range  of  different  solubilities 
including  infinite  solubility  (H'  =  0). 

EXPERIMENT  1 

Experiment  1  was  conducted  on  January  11,  1973.  At  0100  an  occluded 
frontal  system  had  passed  Quillayute  moving  eastward.  At  1000  a  very  deep 
cold,  low  pressure  area  was  located  630  nautical  miles  west-northwest  of 
the  test  site,  oscillating  and  moving  slowly  northward.  The  low  sloped 
west-northwest  aloft,  causing  a  strong  southwesterly  flow  of  cold,  moist, 
unstable  air  with  a  series  of  low-amplitude  short  waves  moving  rapidly 
around  the  low  system  and  onto  the  Washington  and  Oregon  coasts.  This 
brought  shower  activity  and  orographic  precipitation  to  the  test  area.  A 
110-knot  jet  stream  oriented  southwest  to  northeast  over  Hoquiam, 

Washington  was  moving  southward  keeping  the  primary  storm  track  south  of 
the  test  area.  The  freezing  level  was  at  approximately  5200  feet. 


0 


All  precipitation  and  vapor  samplers  were  in  operation  by  1130  PST. 

©  Release  of  ethylacetoacetate  commenced  at  1150,  was  interrupted  momentar¬ 

ily  from  1158  to  1159  when  the  rain  stopped  and  terminated  at  1209.  The 
average  rainfall  rate  was  0.20  cm  h  ;-1  during  the  test  period.  Winds  were 
generally  out  of  the  south-southeast  (157° )  at  about  9  mph  with  occasional 
©  gusts  up  to  20  mph.  A  total  release  of  9.8  liters  was  accomplished  over 

the  18-minute  experiment. 

© 
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TABLE  6. 

OBSERVED  TRACER  CONCENTRATION  VERSUS 
CONCENTRATIONS  PREDICTED  FOR  VARIOUS 
SOLUBILITIES  -  EXPERIMENT  1  (EAA) . 


0 

Station 

^obs . 

CH'»0 

CH'-0.25 

CH’-4 

A  3 

<0.10  ppm 

0.04  ppm 

0.04  ppm 

0.02  ppm 

4 

0.16 

0.19 

0.18 

0.10 

5 

0.50 

0.35 

0.71 

0.34 

© 

6 

1.1 

2.61 

2.41 

0.94 

7 

1.5 

9.10 

7.96 

1.96 

8 

1.6 

36.2 

27.7 

2.97 

9 

1.5 

16.6 

13.9 

2.63 

• 

10 

1.2 

4.98 

4.48 

1.44 

11 

1.0 

0.79 

1.12 

0.52 

12 

0.32 

0.26 

0.25 

0.14 

• 

13 

u.  10 

0.05 

0.05 

0.02 

B  8 

— 

0.06 

0.05 

0.02 

9 

— 

0.24 

0.21 

0.07 

10 

— 

0.80 

0.67 

0.20 

• 

11 

1.0 

2.11 

1.68 

0.43 

12 

0.70 

4.32 

3.18 

0.71 

13 

1.2 

6.39 

4.43 

0.91 

14 

0.78 

5.90 

4.21 

0.87 

• 

15 

0.73 

3.29 

2.54 

0.58 

16 

0.33 

1.48 

0.97 

0.27 

17 

0.35 

0.31 

0.27 

0.09 

18 

0.12 

0.07 

0.06 

0.02 

I 


» 


-30- 


EXPERIMENT  2 


Experiment  2  was  performed  on  the  afternoon  of  January  11,  1973.  The 
synoptic  pattern  and  upper  air  flow  was  therefore  the  same  as  described 
for  Experiment  1.  A  weak  frontal  system  was  developing  approximately  240 
nautical  miles  off  the  Washington  coast  and  most  stations  in  western 
Washington  were  reporting  rain  or  rain  shower  activity.  The  freezing 
level  at  the  time  of  release  had  fallen  to  4200  feet. 

Anemometers  at  the  test  site  showed  light  southeasterly  winds  at  the 
release  height  and  ground  level  winds  with  a  more  easterly  orientation. 

In  response  to  this  observation  three  additional  rain  samplers  were  added 
to  the  western  edge  of  both  arc  A  and  b  to  insure  plume  enclosure.  Release 
of  ethylacetoacetate  was  initiated  at  1519  and  terminated  at  1543  after  a 
discharge  of  11.4  liters  of  the  tracer.  The  rainfall  rate  averaged  approx¬ 
imately  0.59  cm  hour  Winds  were  generally  out  of  the  southeast  (146°) 
at  7  mph  and  gusting  to  a  maximum  of  22  mph. 
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Figure  11.  MEASURED  CONCENTRATIONS  OF  EAA  IN  RAINWATER,  RUN  2. 
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TABLE  /. 

OBSERVED  TRACER  CONCENTRATION  VERSUS 
CONCENTRATIONS  PREDICTED  FOR  VARIOUS 
SOLUBILITIES  -  EXPERIMENT  2  (EAA). 


at  ion 

C  . 
obs . 

CH'«0 

CH'*0. 25 

CH ' 1 

1 

<0.10  ppm 

0.29  ppm 

0.27  ppm 

0.14 

2 

<0.10 

0.94 

0.89 

0.43 

3 

0.80 

2.6 

2 . 38 

1.03 

4 

1.0 

5.8 

5.23 

1.97 

5 

1.5 

12.0 

9.99 

2.94 

6 

2.0 

50.3 

28.7 

3.39 

7 

1.6 

12.0 

9.99 

2.94 

8 

0.77 

5.8 

5.23 

1.97 

9 

0.20 

2.6 

2.38 

1.03 

10 

0.11 

0.94 

0.89 

0.43 

11 

<0.10 

0.29 

0.27 

0.1- 

12 

<0.10 

0.07 

13 

<0.10 

0.02 

5 

<0.10 

0.01 

6 

0.26 

0.05 

7 

0.16 

0.19 

0.16 

0.06 

8 

0.31 

0.62 

0.49 

0.17 

9 

0.36 

1.7 

1.23 

0.38 

10 

0.31 

3.7 

2.44 

0.68 

11 

0.27 

6.2 

3.73 

0.97 

12 

0.37 

7.5 

4.31 

1.09 

13 

0.28 

6.2 

3.73 

0.97 

14 

0.12 

3.7 

2.44 

0.68 

15 

<0.10 

1.7 

1.23 

0.38 

16 

'0.10 

0.62 

0.49 

0.17 

17 

<0.10 

0.19 

0.16 

0.06 

18 

<0.10 

0.05 

* 
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EXPERIMENT  3 


This  experiment  was  conducted  on  January  12,  1973.  At  1400  PST  a  primary 
low  pressure  system  had  moved  northward  into,  the  southwestern  Gulf  of 
Alaska.  A  broad  surface  trough  line  extended  southward  from  the  low  with 
a  developing  wave  in  this  trough  line  located  approximately  480  nautical 
miles  southwest  of  the  test  site  and  moving  to  the  northeast.  A  weak  warm 
front  from  this  wave  intersected  the  Oregon  coast  between  Newport  and 
Astoria  moving  northward  at  27  knots  and  causing  light  overrunning  precip¬ 
itation  to  the  north.  The  upper-air  low  pressure  area  had  moved  eastward 
to  about  700  nautical  miles  off  the  coast  of  Vancouver  Island  causing  a 
lower- atmosphere  flow  of  cold,  moist,  unstable  air  from  the  southwest 
(230°).  Light  shower  activity  and  orographic  precipitation  combined  with 
the  overrunning  to  give  continuous  precipitation  of  variable  intensity. 

The  early  morning  freezing  level  was  at  6300  feet  and  dropped  to  about 
2200  feet  by  1600  PST.  The  main  storm  track  was  kept  to  the  south  of  the 
test  site  by  a  90-knot  jet  stream  oriented  at  260°  and  passing  over 
Newport,  Oregon. 

Prior  to  release  of  diethylamine,  winds  at  the  100  foot  level  were  observed 
to  be  out  of  the  south  (175°)  at  10  mph.  Surface  winds  were  light  and 
slightly  variable.  In  order  to  be  assured  of  total  plume  enclosure,  two 
rain  samplers  were  added  to  the  eastern  edge  of  arc  A  and  three  to  the 
edge  of  arc  B.  Tracer  release  was  initiated  at  1440  and  terminated  at 
1520.  Rainfall  was  variable  throughout  the  test  period  becoming  very  light 
at  1455  but  increasing  significantly  by  1512.  The  average  rainfall  rate 
was  0.25  cm  hour  Winds  were  generally  from  the  south-southeast  (168°) 
at  an  average  speed  of  8  mph  and  gusting  to  26  mph.  A  total  volume  of 
16.6  liters  of  diethylamine  was  dispersed  during  the  40  minute  test 
period. 
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TABLE  8. 

0  OBSERVED  TRACER  CONCENTRATION  VERSUS 

CONCENTRATIONS  PREDICTED  FOR  VARIOUS 
SOLUBILITIES  -  EXPERIMENT  3  (DEA) . 


Station 

^obs . 

CH'-0 

CH’-11 

CH'-21 

0 

A  4 

0.6  ppm 

5 

0.7 

6 

2.1 

0.005  ppm 

0 

7 

1.4 

0.1  ppm 

0.042 

0.013  ppm 

8 

1.3 

1.5 

0.225 

0.084 

9 

2.1 

6.9 

0.628 

0.221 

10 

2.5 

2.2 

0.898 

0.304 

0 

11 

1.0 

6.9 

0.628 

0.221 

12 

0.6 

1.5 

0.225 

0.084 

13 

<0.1 

0.1 

0.042 

0.018 

14 

<0.1 

0.005 

0 

B  11 

0.9 

0.01 

0.002 

0.001 

12 

1.7 

0.12 

0.021 

0.012 

13 

2.2 

0.81 

0.112 

0.058 

14 

2.7 

2.9 

0.304 

0.154 

0 

15 

3.0 

4.8 

0.425 

0.212 

16 

2.1 

2.9 

0.304 

0.154 

17 

2.8 

0.81 

0.112 

0.058 

18 

1.1 

0.12 

0.021 

0.012 

o 

19 

1.2 

0.01 

0.002 

0.001 

0 


a 


u 


EXPERIMENT  4 


The  second  experiment  incorporating  the  use  of  diethylamine  was  conducted 
on  January  16,  1973.  A  deep  low  pressure  system  was  centered  400  nautical 
miles  northwest  of  the  test  site  moving  eastward  at  about  25  knots. 

A  trough  extended  south-southeast  from  the  low  to  200  miles  off  the 
Washington  coast  and  a  small  frontal  wave  over  northcentral  Oregon  was 
filling  and  moving  northeastward.  The  flow  of  cold,  moist,  unstable  air 
aloft  was  southwesterly  bringing  a  series  of  fast  moving,  low- amplitude, 
short  waves  to  the  Washington  coast  and  causing  showers.  A  130-knot  jet 
stream  from  the  southwest  was  located  between  Salem  and  Eugene,  Oregon. 

The  freezing  level  was  at  3800  feet. 

Release  of  diethylamine  was  initiated  at  1129  and  terminated  at  1159.  The 
30-minute  test  was  characterized  by  southerly  winds  (172°)  at  about  16  mph 
and  gusting  to  25  mph.  The  rainfall  rate  was  fairly  uniform  averaging 
0.25  cm  hour  *.  Seventeen  liters  of  the  tracer  were  vaporized  during 
the  test. 
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MEASURED  CONCENTRATIONS  OF  DEA  IN  RAINWATER,  RUN  4. 


TABLE  9 


OBSERVED  TRACER  CONCENTRATION  VERSUS 
CONCENTRATIONS  PREDICTED  FOR  VARIOUS 
SOLUBILITIES  -  EXPERIMENT  4  (DEA). 


Cobs .  CH'-0  CH’-11 


Station 


EXPERIMENT  5 


This  experiment  was  also  conducted  on  January  16,  1973.  The  surface  low 
pressure  system  noted  before  (Experiment  4)  had  moved  to  150  nautical  miles 
off  Vancouver  Island.  The  surface  trough  line  had  filled  and  broadened 
and  the  frontal  wave  locnr.ed  over  north  central  Oregon  was  dissipating 
over  extreme  southeast  Washington.  Other  features  remained  virtually 
unchanged. 

The  generation  of  diethylamine  vapor  was  begun  at  1309.  At  this  time 
heavy  rain  was  falling  which  for  the  first  3-4  minutes  was  mixed  with 
sleet  or  hail.  Winds  at  the  release  height  were  southerly  while  those 
at  the  surface  were  from  the  south-southwest  (,v210o).  The  rainfall  rate 
remained  fairly  constant  throughout  most  of  the  release  period  (1309-1336) 
decreasing  to  a  light  shower  at  1332.  The  average  rainfall  rate  over  the 
course  of  the  experiment  was  1.1  cm  hour  ^  and  winds  were  generally  ob¬ 
served  to  be  from  the  south  (173°)  at  about  16  mph.  The  total  volume  of 
tracer  released  during  this  period  was  7.7  liters. 
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Figure  14.  MEASURED  CONCENTRATIONS  OF  DEA  IN  RAINWATER,  RUN  5. 


'WSTW' 


*  1  »JPJ|  11  r 


TABLE  10. 

OBSERVED  TRACER  CONCENTRATION  VERSUS 
CONCENTRATIONS  PREDICTED  FOR  VARIOUS 


SOLUBILITIES 

-  EXPERIMENT 

5  (DEA). 

Station 

^obs . 

CH'-0 

CH'-11 

CH'-21 

A  8 

0.17  ppm 

9 

0.47 

10 

0.37 

0.04  ppm 

0.02  ppm 

0.01  ppm 

11 

0.43 

0.80 

0.24 

0.10 

12 

0.57 

10.7 

0.70 

0.25 

13 

0.32 

0.80 

0.24 

0.10 

14 

0.48 

0.04 

0.02 

0.01 

15 

0.17 

B  12 

0.28 

13 

0.41 

14 

0.47 

15 

0.47 

0.001 

0.001 

16 

0.47 

0.08 

0.01 

0.019 

17 

— 

0.76 

0.054 

0.112 

18 

0.63 

5.03 

0.097 

0.207 

19 

0.56 

0.76 

0.054 

0.112 

20 

0.43 

0.08 

0.010 

0.019 

21 

0.43 

0.001 

0.001 

22 

0.42 
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CHAPTER  VI. 

DISCUSSION  OF  RESULTS 

As  indicated  by  the  data  in  Tables  6  through  10  ,  the  EPAEC  model 
predicts,  for  the  best  available  estimates  of  solubility,  results  higher 
than  those  measured  for  the  case  of  ethylacetoacetate  and  lower  for  that 
of  diethylamine.  Furthermore,  the  magnitudes  of  the  deviations  are  such 
that  the  EPAEC  model,  with  the  input  parameters  employed,  provides  a  rather 
unsatisfactory  quantitative  description  of  observed  behavior.  We  attrib¬ 
ute  this  lack  of  agreement  primarily  to  the  previously  noted  uncertainties 
in  the  solubility  parameters;  the  pr<_  ent  chapter  presents  an  analysis  of 
the  various  errors  associated  with  this  experiment  in  support  of  this 
conclusion. 

The  possible  sources  of  disagreement  between  predicted  and  experimental 
washout  concentrations  can  be  divided  into  three  categories.  These  are 

•  Experimental  errors 

•  Inadequacies  of  the  washout  model 

•  Inappropriate  values  of  model  input  parameters 

Several  of  these  features  can  be  analyzed  by  noting  behavior  observed 
during  the  EPA's  SO2  studies  conducted  previously  at  the  same  location 
using  the  same  techniques.  These  categories  will  be  addressed  individually 
in  the  following  text. 

EXPERIMENTAL  ERRORS 

The  primary  sources  of  possible  experimental  error  are  itemized  as  follows: 

•  Rain  sampling  and  analysis, 

•  Wind  measurement, 

•  Raindrop  size  measurement, 

•  Release  rate  measurement, 

•  Release  state  of  material 

and  •  Desorption  and  dry  deposition  of  material  on  collectors 

The  errors  associated  with  these  measurements  were  all  small  compared  to 

the  deviations  exhibited  in  the  previous  chapter.  As  indicated  previously 
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the  sampling  and  chemical  analysis  Involves  estimated  errors  of  ^30 
percent  for  dlethylamine  and  <25  percent  for  ethylacetoacetate.  Wind 
measurement  using  the  precisely  leveled  and  calibrated  Gill  Anemometer  was 
not  expected  to  exceed  5 Z  error,  even  under  adverse  conditions;  measure¬ 
ments  of  organic  release  rates  were  valid  to  within  one  or  two  percent, 
Measure',  raindrop  size  spectra  were  open  to  some  question  because  of 
fluctuations  in  rain  characteristics  during  the  course  of  the  experiments. 
Attempts  to  compensate  this  effect  were  made  by  choosing  representative 
spectra  obtained  from  several  measurements  made  during  ei.  h  run.  Although 
there  was  a  fair  degree  of  uncertainty  in  the  "representative"  spectra 
chosen  in  this  manner,  this  was  rendered  of  lesser  importance  owing  to  the 
rather  weak  dependence  of  washout  on  drop  size  expected  under  these  test 
conditions.  Associated  deviations  between  experimental  and  predicted  con¬ 
centrations  are  expected  to  be  less  than  ten  percent. 

Finally,  the  release  state  of  the  organic  material  is  of  some  concern  in 
that  the  plume  was  assumed  (by  the  EPAEC  model)  to  be  in  the  vapor  state 
at  all  times.  This  condition  was  satisfied  well  in  most  cases  by  the 
vapor  generator.  For  ethylacetoacetate,  however,  a  visible  plume  of  fine 
(probably  submicron)  droplets  was  observed  to  exist  for  a  short  distance 
downwind  of  the  generator.  Because  of  the  short  length  of  the  visible 
plume  prior  to  complete  droplet  evaporation  we  estimate  thnt  any  resulting 
deviation  was  less  than  ten  percent  for  the  inner  arc  and  negligible  for 
the  outer  arc. 

An  aspect  of  sampling  not  included  with  the  above  estimates  of  sampling 
and  analysis  errors  involves  the  possibility  of  tracer  desorption  from  the 
rainwater  during  its  time  of  residence  on  the  sampler  funnel  surface. 

This  effect  was  studied  during  the  previous  S02  tests,  and  it  was 
concluded  that  associated  errors  decrease  with  increasing  rain  rate, 
and  the  potential  error  resulting  from  desorption  increases  with 
decreasing  solubility.  For  SO2  (solubility  comparable  to  that  of 
ethylacetoacetate)  at  Quillayute,  this  error  is  within  twenty  percent 
under  the  most  adverse  conditions.  For  lower  solubility  substances 
such  as  dlethylamine,  however,  the  associated  error  is  expected  to  be 
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larger.  Rather  paradoxically,  however,  the  observed  deviations  between 
experiment  and  theory  occur  in  a  direction  opposite  to  that  expected  from 
the  desorption  effect.  This  being  the  case,  one  must  conclude  that  desorp¬ 
tion  effects  were  not  the  primary  source  of  disagreement. 

Two  aspects  of  the  chemical  analysis  of  field  samples  are  also  worthy  of 
mention  at  this  point.  The  apparent  tracer  concentration  of  the  background 
samplers  has  been  alluded  to  previously  in  this  report.  It  is  the  feeling 
of  the  authors  that  particulate  matter  within  the  background  rain  samples 
is  responsible  for  transmittance  values  which  are  generally  less  than  100%. 
Thus  light  scattering  rather  than  absorption  may  account  for  a  fraction  of 
the  decrease  in  transmittance  registered  for  the  field  samples.  Although 
experimental  support  for  this  conclusion  is  lacking,  the  tracer  concentra¬ 
tions  suggested  by  the  analysis  of  background  samples  is  considerably  less 
than  those  observed  for  samples  taken  from  beneath  the  plume.  Error  in¬ 
troduced  by  this  background  material  is  thus  of  such  a  magnitude  as  to  offer 
only  a  minor  influence  of  the  discrepancy  between  the  observed  and  predicted 
sample  concentrations.  Examination  of  the  data  contained  in  Tables  13,  15, 
16,  17,  and  19  (Appendix  D)  confirms  that  the  apparent  tracer  concentration 
of  the  background  samples  rarely  exceeds  20%  of  the  level  found  in  samples 
from  beneath  the  plume. 

In  addition  to  vapor  capture  by  raindrops  passing  through  the  plume  (wash¬ 
out)  dry  deposition  of  the  organic  vapor  at  ground  level  may  offer  a 
-.econd  mechanism  for  tracer  accumulation.  Some  estimate  of  the  importance 
of  this  process  may  be  ventured  by  an  examination  of  the  results  obtained 
from  t'  "*  vapor  samplers.  For  example,  if  the  concentrations  of  vapor  in 
the  gas  ihase  adjacent  to  the  surface  are  demonstrated  to  be  much  lower 
ti  an  those  aloft,  it  is  permissible  (under  the  present  test  conditions)  to 
neglect  the  contribution  of  dry  deposition  relative  to  that  of  scavenging. 
Uaca  acquired  from  the  analysis  of  the  samples  taken  during  Experiment  4 
suggests  that  this  is  indeed  the  case.  These  data  are  tabulated  in  Table  11. 
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TABLE  11. 


RELATIVE  MASS  OF  DIETHYLAMINE  IN  RAIN  AND 
VAPOR  SAMPLES  ACQUIRED  DURING  EXPERIMENT  4. 


Mass  of  DEA 

Mass  of  DEA  in 

Integrated  Vapi 

Station 

in  Rain 

Vapor  Sampler 

Concentration 

A- 7 

21  wg 

1.0  yg 

0.033  yg/1 

8 

22 

1.5 

0.050 

9 

33 

1.7 

0.057 

10 

29 

1.0 

0.033 

11 

28 

<0.8 

<0.027 

12 

20 

<0.8 

<0.027 

13 

19 

<0.8 

<0.027 

14 

11 

<0.8 

<0.027 

B-12 

7.1 

1.3 

0.043 

13 

15 

<0.8 

<0.027 

14 

30 

1.5 

0.050 

15 

44 

2.6 

0.087 

16 

39 

3.6 

0.12 

17 

57 

2.6 

0.087 

18 

58 

2.7 

0.090 

19 

49 

2.4 

0.080 

20 

34 

1.0 

0.033 

In  view  of  the  above  considerations  we  conclude  that  experimental  errors 
are  not  responsible  for  a  significant  portion  of  the  observed  discrepancies 
between  experimental  and  predicted  washout  concentrations.  Subsequent 
sections  of  this  chapter  will  deal  with  additional  aspects  of  this 
disagreement. 
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INADEQUACIES  01  THE  WASHOUT  MODEL 

The  assumptions  employed  during  the  formulation  of  the  EPAEC  model  are 
numerous  and  have  been  described  in  detail  elsewhere  (Dana,  et  al  1973). 
These  Include  the  following: 

•  A  steady  state  exists 

•  The  plume  is  gaseous  and  is  described  by  the  bivariate- 
normal  equation 

•  The  rain  is  composed  of  spherical  droplets,  which  fall 
vertically  with  known  terminal  velocities,  and  move 
laterally  with  the  mean  wind  velocity. 

Our  previous  work  for  the  Environmental  Protection  Agency  on  SO2  washout 
(Dana  1972,  1973)  has  included  a  rather  extensive  analysis  of  these 
assumptions.  Part  of  this  analysis  has  involved  comparisons  of  model- 
predicted  washout  concentrations  with  those  observed  experimentally  under 
a  large  variety  of  circumstances.  In  general  these  comparisons  exhibit 
good  agreement;  washout  distributions  such  as  those  shown  in  Figure  15 
which  pertain  to  an  SO2  experiment  conducted  during  the  EPA  SO2  series, 
are  obtained.  This  evidence  suggests  strongly  that  the  tests  and  basic 
model  described  in  the  present  report,  which  are  in  most  aspects  identical 
to  those  applied  during  the  previous  SO2  study,  should  be  valid  well 
within  the  margins  of  disagreement  exhibited  in  Tables  6  through  10. 

One  can  argue  that  effects  such  as  f luctuating-plume  behavior  and  short¬ 
term  variations  in  rain  rate  should  combine  to  produce  anomalies  not 
considered  here  by  the  model  assessment,  and  thereby  cause  deviations 
between  experiment  and  theory.  While  this  is  certainly  true,  the  fact 
remains  that  markedly  superior  agreement  has  been  obtained  for  the  SO2 
calculations  using  essentially  the  same  modeling  approach  as  that  employed 
here.  In  view  of  this  finding  we  conclude  that  the  primary  source  of 
error  in  the  present  study  does  not  involve  the  basic  structure  of  the 
washout  model. 
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INAPPROPRIATE  VALUES  OF  MODEL  INPUT  PARAMETERS 

The  validity  of  many  of  the  model  Input  parameters  has  been  discussed 
previously  In  the  Experimental  Errors  Section,  above.  Molecular  proper¬ 
ties  (solubility,  dlffuslvlty,  and  viscosity  data),  however,  were  deter¬ 
mined  from  literature  values  or  by  calculation  and  deserve  additional 
consideration  here.  The  kinematic  viscosity  of  air,  used  In  mass-transfer 
calculations,  is  known  with  high  accuracy,  and  therefore  can  be  disregarded 
from  further  consideration  In  this  context.  Dlffuslvlty  estimates,  obtained 
from  collision  integral  theory  as  shown  in  Appendix  B  are  expected  to  be 
in  error  by  no  more  than  ten  percent  (Reid  and  Sherwood  (1958)).  These 
values  are  thus  sufficiently  accurate  to  eliminate  any  concern  for  high 
associated  errors  in  the  washout  calculations. 


As  indicated  previously,  solubility  estimates  were  not  nearly  as  high  in 
quality  as  those  of  the  other  parameters.  These  values,  calculated  as 
shown  In  Appendix  A,  are  open  to  question  on  several  counts.  In  the  case 
of  ethylacetoacetate,  the  solubility  recorded  In  the  Handbook  of  Chemistry 
(Lange,  1956)  Is  not  documented  and  the  conditions  under  which  It  was 
determined  are  thus  unknown.  The  maximum  concentration  of  tracer  within 
the  plume  at  300  feet  is  estimated  by  the  Gaussian  plume  diffusion  model 
to  be  ''-lO  g/liter  for  the  conditions  observed  during  Experiments  1  and 
2.  Since  early  gas  solubility  work  was  generally  directed  to  much  higher 
concentration  levels,  the  handbook  value  is  probably  not  valid  for  the 
test  conditions.  Indeed  the  washout  model  is  only  able  to  predict  rain 
concentrations  approaching  those  observed  If  one  assumes  a  Henry's  Law 
constant,  H',  nearly  twenty  times  larger  than  the  estimated  value  of 

0.25^S£  . 

mole 


Contrary  to  ethylacetoacetate,  a  Henry's  Jaw  constant  has  been  reported 
for  dlethylamlne  although  this  value  Is  also  the  product  of  a  study  at 
rather  high  vapor  concentrations.  There  Is  also  some  question  regarding 
the  vapor  pressures  used  to  adjust  H'  to  the  experimental  conditions.  The 
vapor  pressure  data  of  Stull  (1947)  are  in  conflict  with  that  used  in  the 
calculations  in  Appendix  A.  Stull's  data  would  give  an  H'  value  of 
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8.5  therefore  predict  rein  concentrations  closer  to  those 

observed. 


While  the  ab>  mce  of  high  quality  solubility  data  does  not  permit  an 
assessment  ox  the  EPAEC  model's  ability  to  predict  the  washout  of  spar¬ 
ingly  soluble  organics,  the  data  acquired  from  this  test  series  may  be 
used  to  estimate  the  fractional  washout  per  unit  distance  from  the  source. 
The  mass  of  contaminant  captured  per  kilometer  per  minute,  w.  Is  calculated 
from  rhe  equation 

V4YEC 

w  .  -  -  x  10 


where  I  C*  sum  of  tracer  concentrations  across  the  sampling  arc  In  ppm 

3 

V  ■  average  volume  of  rainwater  samples  in  cm 

2 

A  -  cross-sectional  area  of  funnel  in  cm 
and  t  -  length  of  experiment  in  minutes 

In  order  to  complete  calculations  from  the  data  accumulated  for  each 
sampling  arc,  estimates  of  the  concentrations  of  missing  samples  were  made 
from  the  tracer  loads  observed  in  neighboring  samples.  The  percent  of 
tracer  washed  out  per  kilometer  is  calculated  from  w  and  a  knowledge  of 
the  source  strength,  Q.  A  summary  of  pertinent  field  data  and  the  esti¬ 
mated  fractional  washout  per  kilometer  are  contained  in  Table  12. 

It  la  Interesting  to  note  the  more  efficient  capture  of  diethylamine  by 
rainwater  although  Its  solubility,  as  estimated  from  available  data,  is 
somewhat  lower  than  that  of  ethylacetoacetate.  One  Is  thus  led  to  suspect 
the  efficiency  of  ethylacetoacetate  vapor  generation.  The  low  volatility 
of  this  material  and  the  greater  release  rates  calculated  from  the  field 
notes  lend  support  to  this  contention.  The  lack  of  a  correlation  between 
the  fractional  washout  calculated  for  arcs  A  and  B  is  however  sumrising 
since  one  might  expect  some  capture  of  the  liquid  tracer  by  samplers  In 
Arc  A  under  conditions  of  incomplete  vaporization  and  therefore  a  much 
larger  fractional  washout  from  Arc  A  data  as  opposed  to  that  from  Arc  B. 
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LY  OF  FIELD  DATA  AND  ESTIMATED  FRACTIONAL  WASHOUT  OF  TRACERS 
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CHAPTER  VII. 

CONCLUSIONS 

Rain  scavenging  of  organic  vapors  of  intermediate  solubility  has  been 
investigated.  Comparison  of  rainwater  concentrations  predicted  by  the 
EPAEC  model  with  those  observed  shows  generally  poor  agreement.  This 
comparison  and  an  assessment  of  the  errors  contributed  by  various  input 
parameters  and  the  experimental  design  permit  the  following  conclusions 
to  be  drawn; 

1.  The  assumption  of  infinite  solubility  is  unjustified  for 
the  organic  materials  ethylacetoacetate  and  diethylamine. 
Concentrations  which  are  decidedly  below  those  predicted 

on  the  basis  of  this  assumption  signal  reversible  rainwater- 
organic  vapor  interaction  in  which  pickup  of  the  contaminant 
is  impeded  or  even  reversed  under  some  circumstances. 

2.  The  errors  associated  with  tracer  solubilities  prevent  an 
adequate  appraisal  of  the  predictive  power  of  the  washout 
model.  Solubilities,  which  have  been  estimated  from  the 
available  literature,  must  be  determined  for  the  low  vapor 
concentrations  encountered  in  the  field  experiments  before 
any  conclusions  may  be  reached  regarding  the  adequacy  of 
the  model. 

3.  The  data  derived  from  the  field  study  is  judged  to  be  of 
high  quality.  Errors  contributed  by  the  analytical  tech¬ 
nique,  vapor  generation,  meteorological  data  acquisition 
and  rainwater  collection  system  have  all  been  appraised. 

In  view  of  these  considerations,  it  appears  that  a  thorough 
evaluation  of  the  EPAEC  model  is  contingent  primarily  on 
the  attainment  of  improved  solubility  data. 


ft 


It 
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NOMENCLATURE 


a  drop  radius,  1 

3 

j  concentration  in  liquid  phase  moles/1 

h  plume  release  height,  1 

3 

H'  solubility  parameter,  1  /mole,  cf.  Table  5 

2 

k  mass  transfer  coefficient  moles/1  t 

N  number  of  observations  in  a  wind  sample 

Q  emission  release  rate  (volumetric)  at  source  exit 

3 

conditions  1  /t 

u,v,w  wind  speeds  in  x,y,z  directions  1/t 

vz  raindrop  velocity  in  z-direction  1/t 

x,y,z  distance  coordinates  in  downwind,  crosswind,  and  vertical 

directions,  1 

y  gas-phase  mole  fraction 

C  mass  transfer  parameter  defined  by  equation  (2) 

<J>,0  polar  coordinates  of  wind  direction,  radians 


subscripts 

A  "pollutant  A" 

b  bulk 

x  liquid  phase,  x-direction 

y  gas-phase,  y-direction 

z  vertical  direction 

0,<J>  0,$  directions 


APPENDIX  A 


SOLUBILITY  ESTIMATES 


o 


o 


o 


DIETHYLAMINE 

Solubility  of  diethylamine  was  estimated  from  the  high- concentration 
data  of  Dailey  and  Felsing  (1939),  who  give 


P/c  -  H' 


16.25 


mm  Hg  liter 
mole 


21.38 


(Atm)  cc 
mole 


at  25°  C. 


By 


assuming  applicability  of  the  Clausius-Claperyon  Equation  and  applying 
the  vapor  pressure  data  of  Perry  (1966),  one  can  interpolate  to  estimate 
solubility  at  other  temperatures.  At  5°  C 


21.  38^H~j*  11.28 


ETHYLACETOACETATE 

Lange  (1956) gives  C  *  13  gm  EAA/100  gm  water  at  17°  C. 


a  8.85  x  10 


-4  moles 
cc 


9 


From  Perry  0.963)  the  vapor  pressure  of  EA/  at  17°  P  =  .440  mm  Hg. 

■  5.79  x  1C  ^  Atm 


P 

C 


.65  iAtsJcc 

mole 


at  17°  C 


Using  the  Clausius  Claperyon  Equation  as  before 


H' 


5 


(Atm)cc 

mole 


9 
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APPENDIX  B 


GAS  DIFFUSION  COEFFICIENT  CALCULATIONS 


o 


o 


Diffusion  coefficients  were  estimated  by  applying  collision  integral 
theory  as  described  by  the  Wilke-Lee  modification  of  Hirschfelder's 
relationship  (cf.  Reid  and  Sherwood  (1958)).  The  basic  equation  is 
given  by 

a  T3/2  [(M  +  M  )/M  M  ]1/2 
Day  •  - = -  where 


P  °12  nd 


0.00214-0.000492 


/M1  +  M2\ 

\M1M2  / 


1/2 


Ml*  M2 
T 
P 
o 


12 


and  0 . 


molecular  weights  of  air  and  vapor 

absolute  temperature 
pressure 

intermolecular  force  constant 
collision  integral. 


was  obtained  from  tabulated  values  upon  applying  the  relationships 


12 


l/2(o1  +  o2) 


1.18  $ 


1/3 


3.617  A 


(For  Air) 


C12 

k  “  k 


o 


~  -  97°  K 


(For  Air) 


D 


~  =  1.21  T. 
k  b 
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where  $  is  the  molar  volume  of  material  at  its  normal  boiling  point  T^. 

Tfc  -  4j3°  K  for  EAA. 

T.  -  329°  K  for  DEA. 

D 

Molar  volumes  were  estimated  using  group  contributions  (Reid  and  Sherwood 
(1958)).  The  result  is 


$  -  151.6  ==-  for  EAA. 
mole 


$  ■  82.3  cm 

mole 


for  DEA. 


Combining  the  above  and  applying  the  tables  gives 


°12  " 

4.96 

A 

for 

EAA 

0 

°12  ' 

4.39 

O 

A 

for 

DEA 

e12/k  " 

230° 

K 

for 

EAA 

e12/k  " 

197° 

K 

for 

DEA 

o 

nd* 

1.29 

for 

EAA 

V 

1.22 

for 

9 

DEA 

10°  C 
10°  C 


.0741  ==-  for  EAA 
sec 

.0956  cm^  for  DEA. 
sec 


* 


APPENDIX  C 


A  tipping  bracket  rain  gauge  was  among  the  instruments  employed  to  furnish 
a  description  of  the  test  conditions.  This  unit  was  modified  to  record 
rainfall  in  0.06  mm  increments.  Recordings  made  during  each  of  the  tracer 
releases  are  reproduced  on  the  following  pages. 
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Figure  IS.  RAIN  GAUGE  RECORD  FOR  EXPERIMENT  1  ETHYLACETOACETATE 

RELEASE  -  1150  TO  1158  AND  1159  TO  1209,  JANUARY  11,  1973. 
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APPENDIX  D 


The  following  pages  contain  tabulations  of  the  tracer  concentrations  In 
rain  samples  collected  during  this  test  series.  The  Impinger  solutions 
taken  from  the  vapor  samplers  operated  during  Experiments  1  and  4  have 
been  analyzed  and  are  also  Included.  These  results  gave  an  Indication  of 
the  significance  of  dry  deposition  and  further  analysis  of  impinger  solu¬ 
tions  was  not  undertaken. 
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TABLE  13. 

EXPERIMENT  1. 

Tracer-Ethylacetoacetate 
Time  of  Release  -  18  minutes 
Source  Strength  -  4.6  x  10~2  mole  sec.“l 
Rainfall  Rate 
0.20  cm  hr“l 


Station  Sample  Vol 


EAA  Cone.  EAA  Cone.  -  Less  Background 


Q 


TABLE  13.  (Continued) 

EXPERIMENT  1. 

Tracer-Ethylacetoacetate 
Time  of  Release  -  18  minutes 
Source  Strength  -  A. 6  x  10-2  mole  sec.-3 
Rainfall  Rate 
0.20  cm  hr-3 


o 


0 

Station 

B-6 

O  B-7 

B-8 

B-9 

B-10 

O  B-ll 

B-12 

B-13 

P-14 

O  B-15 

B-16 

B-17 

B-18 

Q  B-19 

B-20 

B-21 

B-22 

O  8-23 

B  -24 

B-25 

B-26 

3  B-27 

B-28 

B-29 


Sample  Vol. 

EAA  Cone. 

41 

0.24  ppm 

40 

J.  76 

Lost 

— 

Lost 

— 

Lost 

— 

37 

1.1 

41 

0.80 

45 

1.3 

44 

0.88 

43 

0.83 

36 

0.43 

42 

0.45 

5 

0.22 

46 

0.28 

46 

0.18 

43 

0.17 

50 

0.18 

46 

0.19 

46 

0.15 

46 

0.22 

57 

0.08 

53 

0.07 

59 

0.24 

53 

0.15 

EAA  Cone.  -  Less  Background 

0.1A  ppm 
0.66 


1.0 
0.70 
1.2 
0.78 
0.73 
0.33 
0.35 
0.12 
0.18 
<0.1 
<0.1 
<0.1 
<0.1 
<0.1 
0.12 
<0.1 
<0.1 
0.1A 
-0.  1 


TABLE  14 


o 


o 


o 


0 


0 


0 


0 


ft 


ft 


Station 

EXPERIMENT  1. 

Impinger  Solutions 

Mass  of  Ethvlscetoacetate 

BKG 

0.8 

A- 3 

0.8 

4 

0.0 

5 

1.6 

6 

1.2 

7 

0.8 

8 

0.8 

9 

1.9 

10 

0.4 

11 

0.4 

12 

0.2 

13 

0.0 

14 

0.0 

B-10 

2.6 

11 

1.3 

12 

1.6 

14 

2.6 

15 

2.1 

16 

2.1 

17 

1.0 

18 

2.4 

19 

0.8 

20 

1.6 
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TABLE  15. 

EXPERIMENT  2. 

Tracer-E thy lace toace tat e 
Time  of  Release  -  24  minutes 


Source  Strength  -  6.3  x  10“2  mole  sec. 


0 

Rainfall  Rate 
0.59  cm  hr~l 

Sampling  Station 

Sample 

Vol.  EAA  Cone. 

EAA  Cone.  -  Less  Background 

BKG 

105 

ml  0.08  ppm 

0 

AC-3) 

~100 

0.08 

<0.1  ppm 

A(-2) 

87 

0.17 

<0.1 

A(-l) 

~100 

0.08 

<0.1 

A-l 

141 

0.10 

<0.1 

0 

A- 2 

127 

0.05 

<0.1 

A- 3 

137 

0.88 

0.80 

A- 4 

127 

1.1 

1.0 

A- 5 

138 

1.6 

1.5 

0 

A- 6 

132 

2.1 

2.0 

A- 7 

135 

1.7 

1.6 

A- 8 

133 

0.85 

0.77 

A- 9 

133 

0.28 

0.20 

0 

A- 10 

126 

0.19 

0.11 

A- 11 

133 

0.14 

<0.1 

A- 12 

128 

0.07 

<0.1 

A-l  3 

~100 

0.13 

<0.1 

0 

A-14 

~  70 

0.10 

<0.1 

A-15 

121 

0.06 

<0.1 

A- 16 

123 

0.06 

<0.1 

A- 17 

134 

0.05 

<0.1 

• 

B(-3) 

85 

0.18 

<0.1 

B(-2) 

130 

0.07 

<0.1 

B  (-1) 

81 

0.18 

<0.1 

» 

1 

B-l 

127 

0.06 

<0.1 

B-2 

139 

0.15 

<0.1 
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TABLE  15.  (continued) 


EXPERIMENT  2. 

Tracer-Ethylacetoacetate 
Time  of  Release  -  24  minutes 
Source  Strength  -  6.3  x  10-2  mole  sec.-* 


Sampling  Station 

Rainfall  Rate 
0.59  cm  hr“l 
Sample  Vol.  EAA  Cone. 

E/A  Cone.  -  Less  Background 

B-3 

131 

0.06  ppm 

<0.1  ppm 

B-4 

138 

0.16 

<0.1 

B-5 

130 

0.16 

<0.1 

B-6 

142 

0.34 

0.26 

B-7 

126 

0.24 

0.16 

B-8 

127 

0.39 

0.31 

B-9 

137 

0.44 

0.36 

B-10 

139 

0.39 

0.31 

B-ll 

126 

0.35 

0.27 

B-12 

140 

0.45 

0.37 

B-13 

133 

0.36 

0.28 

B-14 

114 

0.20 

0.12 

B-15 

129 

0.13 

<0.1 

B-16 

125 

0.04 

<0.1 

B-17 

117 

0.10 

<0.1 

B-18 

108 

0.06 

<0.1 

3-19 

133 

0.04 

<0.1 

B-20 

128 

0.04 

<0.1 

B-21 

111 

0.06 

<0.1 

B-22 

113 

0.06 

<0.1 

B-23 

128 

0.05 

<0.1 

B-24 

130 

0.05 

<0.1 

B-25 

125 

0.05 

<0.1 

B-26 

114 

0.07 

<0.1 

B-27 

— 

— 

B-28 

133 

0.04 

<0.1 

B-29 

130 

0.04 

<0.1 
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TABLE  16. 
EXPERIMENT  3. 


Tracer-Diethyl amine 
Time  of  Release  -  40  minutes 
Source  Strength  -  1.7  x  10~1  mole  sec.-! 
Rainfall  Rate 
0.25  cm  hr-^ 
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Less  Background 


Sampling  Station 

Sample  Vol. 

DEA  Cone. 

DEA  Cor 

• 

BKG 

67  ml 

0.0  ppm 

BKG 

36 

0.3 

A-l 

43 

0.2 

<0.1 

A- 2 

43 

0.3 

0.1 

• 

A- 3 

43 

0.4 

0.2 

A-4 

43 

0.8 

0.5 

A- 5 

39 

0.9 

0.7 

A- 6 

38 

2.3 

2.1 

• 

A- 7 

42 

1.6 

1.4 

A- 8 

41 

1.5 

1.3 

A- 9 

43 

2.3 

2.1 

A-10 

45 

2.7 

2.5 

• 

A-ll 

39 

1.2 

1.0 

A- 12 

44 

0.8 

0.6 

A-l  3 

38 

0.2 

<0.1 

A-14 

39 

0.1 

<0.1 

• 

A-15 

34 

0.9 

0.7 

A-16 

40 

1.1 

0.9 

A-l  7 

41 

0.2 

<0.1 

A- 18 

44 

0.4 

0.2 

• 

A- 19 

30 

1.9 

1.7 

A-20 

~40 

0.3 

0.1 

B-l 

41 

0.3 

0.1 

• 

B-2 

41 

0.2 

<0.1 

i  nr  -r • 


TABLE  16.  (continued) 
EXPERIMENT  3. 


Tracer-Diethylamine 

Time  of  Release  -  40  minute 
Source  Strength  -  1.7  x  10“  1  mole 

3 

sec.-l 

Sampling  Station 

Rainfall 
0.25  cm 

Sample  Vol. 

Rate 

hr"l 

DEA  Cone. 

Less  Background 
DEA  Cone. 

B-3 

39 

0.1  ppm 

<0.1  ppm 

B-4 

44 

0.3 

0.1 

B-5 

42 

0.2 

<0.1 

B-6 

41 

1.2 

1.0 

B-7 

44 

0.5 

0.3 

B-8 

37 

0.6 

0.4 

B-9 

42 

0.5 

0.3 

B-10 

39 

1.1 

0.9 

B-ll 

41 

1.1 

0.9 

B-12 

40 

1.9 

1.7 

B-13 

40 

2.4 

2.2 

K-14 

36 

2.9 

2.7 

B-15 

34 

3.2 

3.0 

B-16 

39 

2.3 

2.1 

B-17 

38 

3.0 

2.8 

B-18 

40 

1.3 

1.1 

B-19 

43 

1.4 

1.2 

B-20 

43 

0.6 

0.4 

B-21 

37 

0.3 

0.1 

B-22 

43 

0.3 

0.1 

B-23 

38 

0.4 

0.2 

B-24 

45 

0.4 

0.2 

B-25 

41 

0.4 

0.2 

B-26 

42 

0.1 

<0.1 

B-27 

39 

0.1 

<0.1 

B-28 

42 

0.2 

<0.1 
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TABLE  16.  (continued) 


EXPERIMENT  3. 

Tracer-Diethylamine 
Time  of  Release  -  40  minutes 
Source  Strength  -  1.7  x  10“  ^  mole  sec.'^ 

Rainfall  Rate 
0.25  cm  hr-* 

Less  Background 


Sampling  Station 

Sample  Vol. 

DEA  Cone. 

DEA  Cone. 

B-29 

34 

0.2  ppm 

<0 . 1  ppm 

B-30 

37 

0.2 

<0.1 

B-31 

34 

0.5 

0.3 

B-32 

32 

0.0 

<0.1 

* 


t 
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TABLE  17 


O 


o 


o 


c 


a 


T, 


EXPERIMENT  A. 

Trace r-Die thy lamine 
Time  of  Release  -  30  minutes 
Source  Strength  -  3.05  x  10"!  mole  see.-! 

Rainfall  Rate 

0.25  cm  hr-*  .  _  .  . 

-Less  Background 

Sampling  Station  Sample  Vol.  DEA  Cone.  DEA  Cone. 


29  ml 


0.29  ppm 


A- 1 

2A 

0.29 

<0.10 

A-2 

31 

0.25 

<n.io 

A- 3 

32 

0.17 

<0.10 

A- A 

23 

0.07 

<0.10 

A- 5 

17 

0.52 

0.19 

A- 6 

33 

0.52 

0.19 

A  -7 

38 

0.88 

0.55 

A-8 

36 

0.95 

0.62 

A-9 

3A 

1.29 

0.96 

A-10 

22 

1.6A 

1.3 

A- 11 

28 

1. 30 

1.0 

A- 12 

23 

1.18 

0.85 

A  -13 

33 

0.89 

0.56 

A-]  A 

28 

0.72 

0.39 

A-.'5 

30 

0.20 

<0.10 

A- 16 

36 

0.11 

<0.10 

A- 17 

36 

1.07 

0.7A 

A- 18 

30 

1.A9 

1.2 

B-l 

3A 

0.52 

O.io 

B-2 

31 

0.51 

0.18 

B-3 

A1 

0.. 

<0.10 

B-A 

32 

0.13 

<0.10 

B-5 

A0 

0.14 

<0.10 
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TABLE  17.  (continued) 


# 


EXPERIMENT  4. 


T racer- Diethy lamine 
Time  of  Release  -  30  minutes 
Source  Strength  -  3.05  x  10-1  mole  sec. 
%  Rainfall  Rate 

0.25  cm  hr“l 


-Less  Background 


Sampling  Station 

Sample  Vol. 

DEA  Cone. 

DEA  Cone. 

B-6 

36 

0.22  ppm 

<0.10  ppm 

h 

B-7 

38 

0.21 

<0.10 

B-8 

37 

0.22 

<0.10 

B-9 

27 

0.36 

<0.10 

B-10 

32 

0.12 

<0.10 

9 

B-ll 

33 

0.28 

o 

o 

V 

B-12 

34 

0.54 

0.21 

B-13 

35 

0.75 

0.42 

B-14 

23 

1.6 

1.3 

9 

B-15 

34 

1.6 

1.3 

B-16 

23 

2.0 

1.7 

B-17 

26 

2.5 

2.2 

B-18 

34 

2.0 

1.7 

» 

B-19 

35 

1.7 

1.4 

B-20 

19 

2.1 

1.8 

B-21 

28 

1.6 

1.3 

B-22 

32 

1.1 

0.80 

9 

B-2  3 

27 

0.78 

0.45 

B-24 

~  35 

0.43 

0.  10 

B-25 

35 

0.72 

0.  39 

B-26 

~  35 

0.25 

<0.  10 

• 

B-27 

32 

0.24 

•'O.IO 

B-28 

34 

0.27 

<0. 10 

B-29 

32 

0.05 

0.10 

» 
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TABLE  18 


EXPERIMENT  4. 

Impinger  Solutions 

Station  Mass  of  Diethylamine 


BKG 

<0.8  pg 

A- 3 

2.0 

A 

<0.8 

5 

<0.8 

6 

1.7 

7 

1.0 

8 

1.5 

9 

1.7 

10 

1.0 

11 

<0.8 

12 

<0.8 

13 

<0.8 

14 

<0.8 

B-9 

<0.8 

10 

<0.8 

11 

<0.8 

12 

1.3 

13 

<0.8 

14 

1.5 

15 

2.6 

16 

3.6 

17 

2.6 

18 

2.7 

19 

2.4 

20 

1.0 

1 

i 
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TABLE  19 


EXPERIMENT  5. 


Tracer-Die thy lamine 
Time  of  Release  -  25  minutes 
Source  Strength  -  2.6  x  10“  1  mole  sec."* 

Rainfall  Rate 
1.1  cm  hr”l 

Less  Background 


Sampling  Station 

Sample  Vol. 

DEA  Cone. 

DEA  Cone. 

BKG 

106  ml 

0.15  ppm 

BKG 

107 

0.07 

A-l 

108 

0.13 

<0.10  ppm 

A- 2 

94 

0.04 

<0.10 

A- 3 

84 

0.13 

<0.10 

A- 4 

63 

0.20 

<0.10 

A- 5 

101 

0.31 

0.20 

A- 6 

109 

0.21 

0.10 

A- 7 

125 

0.25 

0.14 

A- 8 

123 

0.28 

0.17 

A-9 

120 

0.58 

0.47 

A-10 

52 

0.48 

0.37 

A- 11 

95 

0.54 

0.43 

A- 12 

48 

0.68 

0.57 

A-13 

51 

0.43 

0. 32 

A- 14 

85 

0.59 

0.48 

A- 15 

99 

0.28 

0.17 

A-16 

76 

0.  3^ 

0.23 

A- 17 

96 

0.33 

0.22 

A- 18 

79 

0.52 

0.41 

B-l 

80 

0.27 

0.  16 

B-2 

86 

0.19 

<0.10 

B-3 

104 

0.21 

0.10 

B-4 

90 

0.16 

<0.  10 

B-5 

81 

0.44 

0.33 
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TABLE  19.  (continued) 

EXPERIMENT  5. 

Tracer-Die thy lamine 
Time  of  Release  -  25  minutes 
Source  Strength  -  2.6  x  lO-*  mole  soc.-^ 

Rainfall  Rate 
1.1  cm  hr-l 

Less  Background 

Sampling  Station  Sample  Vol.  DEA  Cone.  DEA  Cone. 


0.05  ppm 


<0.10  ppm 


.  4  / 

1  4 

• 

.63 

.56 

.43 

.43 

j 

.42 

.45 

I  1 

.32 

.33 

; 

Jill 


APPENDIX  i: 


DESCRIPTION  OK  COMPUTER  CODE 
FOR  THE  EPAEC  NONLINEAR  NON  FEEDBACK  WASHOUT  MODEL 


The  computer  code  can  be  described,  in  a  somewhat  superficial  sense,  in 
terms  of  the  flow  chart  shown  in  Figure  2.  This  figure  demonstrates  that 
washout  calculations  may  be  performed  using  the  code  simply  by  employing 
a  main  program  that  performs  the  following  functions:  1)  reading  of 
appropriate  input  data;  2)  execution  of  algorithm  by  the  statement  CALL 
MASTER;  3)  printing  of  resulting  computed  values.  MASTER  is  a  master 
coordinating  subroutine  which  employs  all  of  the  program  functions  to 
calculate  washout  as  indicated  in  Figure  2. 

The  utility  of  the  above  arrangement  is  that  it  allows  the  code  to  be 
applied  generally  for  a  variety  of  specific  purposes;  one  simply  writes  a 
main  program  designed  to  fulfill  his  particular  requirements,  and  employs 
the  statement  CALL  MASTER  to  execute  the  basic  algorithms  required.  The 
Centralia  power  plant  computations  described  in  our  previous  report 
(Dana,  et  al  (1973)),  for  example,  were  executed  using  a  main  program 
that  memorized  the  topography  of  the  surrounding  area,  computed  relative 
distances  based  on  the  plume  location,  and  then  calculated  corresponding 
washout  concentrations  by  the  statement  CALL  MASTER. 

A  generalized  main  program  can  be  employed  to  perform  the  above  functions, 
if  desired.  Use  of  such  a  program  will  provide  results  identical  to  those 
computed  using  customized  main  programs,  the  major  disadvantages  being 
probable  increased  inconvenience  in  data  input  and  less  cont'ol  over  output 
formats.  An  example  of  such  a  generalized  program  is  shown  on  the  pages 
immediately  following,  which  define  the  computer  nomenclature.  Named  EPAEC, 
this  program  employs  a  common  statement  to  facilitate  exchange  of  informa¬ 
tion  with  the  subroutine  MASTER.  The  following  EXTERNAL  statement  desig¬ 
nates  the  function  subroutines  V,  HPFTME,  and  YAB  to  be  used  internally; 
this  statement  is  essential  to  the  operation  of  anv  main  program  for  this 
purpose . 
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Reading  of  input  data  proceeds,  governed  by  the  control  variables  J1-J4 
which  allow  the  bypass  of  designated  READ  statements  if  desired.  Finally, 
an  optional,  internally-generated  receptor  spacing  (DYDUM)  is  calculated. 
This  provides  for  performance  of  subsequent  calculations  at  5-degree  in¬ 
creasing  cross-plume  distances. 

Upon  printing  the  input  data,  MASTER  is  called.  Execution  of  this  sub¬ 
routine  pertains  to  a  single  receptor  location.  Upon  completion  of  the 
calculations  controlled  by  MASTER,  appropriate  printing  is  performed,  in¬ 
cluding  that  of  REF,  which  is  the  variable  denoting  height  above  the  source 
passed  by  individual  raindrops  at  x  *  0. 

At  this  point,  depending  upon  the  current  value  of  the  control  variable 
JEND,  the  program  has  two  options.  It  can  either  proceed  to  the  adjoining 
cross-plume  location,  perform  subsequent  computations,  and  integrate  to 
provide  (ultimately)  a  downwind  washout  rate  (WORATE)  or  it  can  read  new 
input  data  and  proceed  to  calculations  for  other  conditions.  In  the  event 
that  cross-plume  integrations  are  performed,  the  program  terminates  when 
relatively  low  concentrations  are  encountered  on  the  edges  of  the  plume, 
and  the  downwind  washout  rate  is  printed. 

The  program  is  terminated  completely  whenever  no  additional  data  cards  are 
found,  or  when  the  control  variable  JEND  is  set  equal  to  1.  An  example 
data  set  for  use  with  EPAEC  has  been  given  previously  in  Chapter  4. 


DESCRIPTION  OF  BASIC  COMPUTATIONAL  ALGORITHM 

Exclusive  of  the  input-output  functions  governed  by  the  main  program,  the 
basic  computational  algorithm  can  be  described  by  the  hierarchy  of  sub¬ 
routines  shown  in  Figure  2.  The  primary  function  of  this  algorithm  is  to 
solve  the  drop-response  equation 


3K 

—  (v  . 
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which  was  given  previously  (Hales,  et  at  (1972)).  The  calculations  per¬ 
formed  by  the  computer  listing  given  here  envision  a  rain-plume  situation 
as  shown  in  Figure  1. 
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In  this  visualization  a  single  raindrop  of  radius  R  falls  through  the  plume 

to  a  receptor  over  a  linear  trajectory  determined  by  the  wind  speed  U  and 

the  terminal  fall  velocity  V.  A  computational  grid  is  set  up  along  the 

trajectory,  and  Equation  El  is  solved  numerically  using  a  Runge-Kutta 

finite-difference  approximation  to  obtain  the  value  of  c.,  (CGRND  in  the 

Ab 

computer  code)  at  the  receptor  location.  This  procedure  is  repeated  for 
all  drop  sizes  in  the  discretized  spectrum  and  the  results  are  averaged  to 
obtain  a  final  mixed-mean  concentration  (CAVG)  for  that  location. 

y..  values  along  the  trajectory  are  famished  by  a  plume  model  in  the  form 
of  a  subroutine  (YAB) .  At  present  this  subroutine  is  written  to  accomodate 
the  Gif ford-Pasquill  Bivariate-Normal  plume  model,  modified  to  account  for 
quasif irst-order  chemical  reaction,  plume  loft,  and  ambient  background. 
Descriptions  of  the  subroutines  shown  in  Figure  2  are  given  individually 
in  the  following  text.  Listings  of  these  subroutines  are  provided  at  the 
end  of  this  Appendix. 

Subroutine  V  provides  the  terminal  fall  velocity  dz/dt,  as  a  function  of 
raindrop  radius  R.  This  internal  function  is  simply  an  empirical  poly¬ 
nomial  fit  to  measured  terminal  velocity  data. 

Subroutine  HPRIME  utilizes  appropriate  computations  thereof,  to  determine 
values  of  the  apparent  Henry' s-law  constant  appropriate  to  the  current 
value  of  either  c^  or  y  depending  on  the  value  of  the  control  variable 
LGOPT.  If  LGOPT  is  set  equal  to  1  in  the  calling  sequence,  calculation  is 
based  on  the  value  of  y  otherwise  the  current  value  of  c^  is  employed. 

Subrout  ne  MTC  computes  the  overall  mass-t rans fer  coefficient  K^.  Tt  begins 
by  calculating  the  gas-phase  coefficient.  Then,  if  the  control  variable 
JOPT  has  been  set  equal  to  1,  (gas-phase  limiting),  the  overall  coefficient 
is  returned  as  the  gas-phase  value.  Otherwise,  a  liquid  coefficient  is 
calculated,  and  the  gas  and  liquid  coefficients  are  combined  to  obtain  the 
corresponding  stagnant-drop  values. 

Subroutine  RUNGE,  in  conjunction  with  subroutine  SOLVE  performs  the  numeri¬ 
cal  integration  of  Equation  El,  SOLVE  supplying  RUNGE  values  of  the 

derivative  dc., /dz  (F  in  the  computer  code)  and  RUNGE  returning  values  of 
Ab 
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.  RUNGE  is  a  fourth-order  Runge-Kutta  algorithm  which  has  been  adapted 
from  a  previous  work.  This  algorithm  is  rather  complex  and  will  not  be 
discussed  in  detail  here,  except  to  state  that  its  expected  errors  are  of 
the  order  of  the  fractional  grid  spacing  to  the  fourth  power.  For  a  com¬ 
plete  discussion  of  this  method  the  reader  is  referred  to  the  work  of 
Carnahan,  et  al  (1969). 

Subroutine  YAB,  as  mentioned  previously,  provides  values  of  the  gas-phase 
mole  fraction  of  pollutant  as  a  function  of  spatial  location  x,  v,  and  z. 

The  Gif ford-Pasquill  Bivariate-Normal  equation,  modified  for  quasifirst- 
order  gas-phase  chemical  reaction,  is  employed  for  this  purpose.  The 
subroutine  computes  values  of  the  dispersion  parameters,  and  proceeds 
directly  with  a  solution  of  the  bivariate-normal  equation  to  provide  the 
return  value  YAB.  H  in  this  subroutine  is  a  virtual  value,  and  depends 
upon  the  value  of  emission  height  supplied  by  the  calling  program. 

Subroutine  SOLVE  performs  the  function  of  establishing  the  computational 
grid  and  implementing  RINGE  to  obtain  the  solutions  to  Equation  El 
The  routine  begins  by  initializing  variables  and  then  testing  for  the 
occurrence  of  equilibrium  scavenging.  This  is  accomplished  by  calculating 
a  virtual  emission  height  (HSTAR)  and  an  appropriate  dispersion  parameter 
(SIGMAZ) .  This  is  followed  by  determination  of  an  effective  Henry's-law 
constant,  wM ''t*  is  employed  in  conjunction  with  the  criterion  (GROUP  <  15) 
to  determine  whether  or  no*-  equilibrium  scavenging  occurs.  If  equilibrium 
conditions  are  indeed  predicted,  the  scheme  bypusses  the  solution  of 
Equation  £1,  and  simply  returns  the  equilibrium  washout  concentration  value. 

In  the  event  that  equilibrium  conditions  are  not  predicted,  the  routine 
initializes  the  concentration  of  pollutant  in  the  drop,  C(l),  to  its 
appropriate  above-plume  value,  and  proceeds  to  establish  a  computation 
grid.  In  performing  this  function  it  first  tests  for  plume  undercut  by 
the  raindrop  (REF  <  0).  In  the  event  that  undercutting  occurs,  the  grid 
network  is  established  by  dividing  the  vertical  distance  between  the 
sampler  and  the  height  at  which  the  drop  crosses  X  =  0  into  thirty  equal lv- 
spaced  measurements. 
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If  plume  undercutting  does  not  occur,  a  "normal"  grid  spacing  is  estab¬ 
lished.  This  is  accomplished  by  finding  an  appropriate  vertical  spread 
parameter  (SIGMA)  for  early  stages  of  the  drop-plume  encounter  and  (rather 
arbitrarily)  beginning  numerical  computations  at  an  elevation  equal  to  the 
effective  release  heights  plus  three  times  the  computed  spread  (HSTAR  +  3 
*  SIGMA).  Grid  spacing  is  set  at  one-sixtieth  of  the  vertical  distance 
between  the  receptor  and  the  point  where  calculations  are  initiated. 

A  final  modification  of  the  comoutation  grid  structure  is  performed  if  the 
raindrop  encounters  a  plume  having  a  low  degree  of  spread  ("compact 
plume").  This  is  done  simply  by  testing  for  whether  the  current  spacing 
is  less  than  one-fourth  the  computed  spread  parameter  SIGMA.  If  not,  a 
top  grid  spacing  (TD2)  is  set  equal  to  SIGMA/4  (SIGMA/8  if  SIGPHI  is  greater 
than  0.5).  If  "compact"  plumes  are  encountered,  this  finer  grid  spacing  is 
employed  for  25  increments,  and  the  original  grid  spacing,  DELTA7 ,  is 
employed  thereafter.  The  choices  of  grid  spacings  described  above  were 
arrived  at  after  experimentation  with  various  arrangements.  This  system 
provides  for  general  stability  and  accuracy  of  the  algorithm,  with  reason¬ 
able  economy  in  execution  time. 

As  described  earlier,  numerical  solution  of  the  object  equation  is  accom¬ 
plished  using  subroutine  RUNGE.  This  subroutine  is  called  repeatedly,  and 
control  is  transferred  between  it  and  the  calling  subroutine  SOLVE,  which 
updates  the  downwind  distance  (x)  and  effective  release  height  HSTAR. 

SOLVE  also  updates  the  derivative  function  F ( 1 )  and  su>plies  this  value 
to  RUNGE,  which  in  turn  provides  calculated  values  of  tne  concentration 
C(l). 

Calculations  continue  until  the  receptor  location  is  encountered.  Then 
the  value  of  the  ground-level  rain  concentration  (COBJ)  is  calculated  and 
returned  with  other  pertinent  variables  in  the  calling  sequence. 

Subroutine  MASTER  coordinates  calculations  for  the  ensemble  of  drops  in 
the  discretized  spectrum,  and  combines  the  resulting  concentration  values 
to  obtain  mixed-mean  levels.  MASTER  simply  calls  subroutine  SOLVE  for 
each  raindrop  size  in  the  speett  .m,  saves  tb.'  individual  concentrations  in 
the  array  CGRND(I),  and  averagei  accruing  tc  the  equation 
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Control  is  then  transferred  to  the  main  program  for  subsequent  printing 
operations . 

MODIFICATION  OF  THE  COMPUTER  CODE 

The  modular  form  of  the  general  computer  code  enables  it  to  be  modified 
easily  for  use  in  other  applications.  Such  modifications  can  be  categorized 
into  two  tvpes,  depending  on  whether  they  are  meant  to  improve  the  accuracy 
of  the  calculations  or  to  adapt  the  algorithm  for  use  with  substances 
other  than  those  originally  intended. 

The  first  type  of  modification — an  incorporation  of  an  improved  plume 
model  for  instance — can  usually  be  accomplished  by  modular  replacement  of 
one  or  more  subroutines  in  a  straight-forward  manner.  The  second  type  of 
modification  usually  can  be  accomplished  easily,  depending  on  the  materials 
of  interest.  If  this  material  is  a  nonreactive  gas,  one  simply  must  re¬ 
place  the  solubility  function  HPRIME  with  one  appropriate  to  the  gas  in 
question.  Other  routines  are  generally  applicable,  and  corresponding 
modifications  are  accomplished  automatically  by  changes  in  the  physical- 
properties  input  data.  Scavenging  of  a  totally  soluble  gas,  for  instance, 
can  be  calculated  simply  by  modifying  the  function  HPRIME  to  return  a  zero 
value  whenever  it  is  called. 


t 
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TABLE  20. 


COMPUTER  NOMENCLATURE 


0 

Symbol 

AA 

Units 

2 

cm 

. 

B 

gm  moles/liter 

o 

BKG 

dimensionless 

C 

gm  moles/ cm  ^ 

o 

CAVG 

gm  moles/cm^ 

CAVGCL 

gm  moles/cm^ 

o 

CCUM 

gm  moles/cm^ 

CDUM 

2 

gm  moles/ (liter)" 

o 

CEQ 

gm  moles/ cm ^ 

CGRND 

gm  moles /cm^ 

o 

CLAST 

gm  moles/cm^ 

COBJ 

gm  moles/cm^ 

3 

CSET 

gm  moles/cm^ 

CTEST 

gm  moles/cm^ 

IT 


Definition 

Summing  variable  for  calculation  of  average 
concentration 

Dummy  variable  in  subroutine  HPRIME  used  for 
storage  of  sum  of  the  equilibrium  constant 
and  the  excess  hydrogen  ion  concentration 

Mixing  ratio  of  pollutant  in  gas  phase  back¬ 
ground  (moles/mole) 

Mixed-mean  concentration  of  pollutant  in  a 
specific  raindrop 

Mixed-mean  concentration  of  pollutant  in  a 
collected  rain  sample 

Mixed-mean  concentration  cf  pollutant  in  rain 
sample  collected  beneam  plu-  ’  centerline 

Cumulative  concentration  used  for  integration 
across  the  plume  to  calculate  washout  rate 

Dummy  variable  used  in  HPRIME 

Concentration  of  pollutant  in  rain  in  equilib¬ 
rium  with  ground-level  gas-phase  concentration 

Mixed-mean  concentration  of  pollutant  in  a 
specific  raindrop  at  receptor 

Dummy  variable  used  for  performing  cross-pium^ 
integration 

Mixed-mean  concentration  of  pollutant  in  a 
specific  raindrop  at  receptor 

Mixed-mean  concentration  of  pollutant  in  a 
specific  raindrop  used  for  calculation  of 
mass-transfer  coefficient  in  subroutine  TRY 

Dummy  variable  used  to  test  for  cross-plume 
integration  termination  conditions;  also  used 
to  establish  compact  plume  characteristics  in 
SOLVE 
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TABLE  20.  (Contd.) 


Symbol 

Units 

Definition 

D  cm 

Raindrop  diameter 

DAX 

2, 

cm  /sec 

Molecular  diffusivity  of  pollutant  in  air 

DAY 

2, 

cm  /sec 

Molecular  diffusivity  of  pollutant  in  water 

DELTAY 

cm 

Cross-plume  spacing  of  calculation  points 

DELTAZ 

cm 

Vertical  grid  spacing  employed  under  non¬ 
compact  plume  conditions 

DYDUM 

cm 

Cross-plume  spacing  of  calculation  points 
generated  by  computer  if  no  value  is  entered 
as  data 

DZ 

cm 

Vertical  grid  spacing 

DZTST 

cm 

Test  variable  for  assessing  compact  plume 
conditions 

EQCON 

gm  moles/liter 

Equilibrium  constant  for  first  dissociation 
of  SO2  (cf.  Equation  (40)) 

F 

-3 

V.  . 

Probability-density  function  for  raindrops  of 
size  R  in  a  distributed  system;  also  denotes 
derivative  used  in, subroutines  SOLVE  and 

RUNGE  (gm-moles/cm  ) 

GR 

dimensionless 

Dummy  variable  used  in  subroutine  SOLVE 

GROUP 

dimensionless 

Dimensionless  group  used  to  evaluate  equilib¬ 
rium  scavenging  conditions  (cf.  Equation  (1)) 

H 

cm 

Effective  emission  release  height 

HANK 

liters/gm  mole 

Henry 's-law  constant  for  undissociated  SO2  in 
water  (cf.  Equation  (39)) 

HEX 

gm  moles/liter 

Hydrogen  ion  in  rain  other  than  that  contrib¬ 
uted  by  dissolved  SO2 

HP  RIME 

cm^/gm  mole 

Effective  Henry's-law  constant  for  total 
dissolved  SO2  in  water 

HSTAR 

cm 

Plume  height,  or  effective  release  height 
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TABLE  20.  (Contd.) 


Symbol  Units 

3 

HTEST  cm  /gm  mole 

I 

ICP 

IGKND 

IS 

IY 

J 

J1 

J2 

J3 

J4 

JEND 

JOPT 

JP 

.GOPT 


Definition 

Effective  Henry's-law  constant  used  for 
evaluation  of  equilibrium  scavenging  conditions 

Index  integer 

Internal  control  variable  in  subroutine  SOLVE 
providing  for  grid  spacing  modifications  in 
the  case  of  a  compact  plume 

Internal  control  variable  in  subroutine  SOLVE 
providing  for  termination  of  the  algorithm  as 
the  raindrop  encounters  the  receptor 

Internal  control  variable  communicating  status 
of  solution  between  subroutines  SOLVE  and 
RUNGE 

Internal  control  variable  providing  a  runaway 
trap  on  cross-plume  integration  sequence 

Index  integer 

Read  control  variable  providing  for  optional 
reading  of  raindrop  size  distribution  data 

Read  control  variable  providing  for  optional 
reading  of  physical  properties  data 

Read  control  variable  providing  for  optional 
reading  of  plume  data 

Read  control  variable  providing  for  optional 
reading  of  grid  data 

Program  termination  control  variable 

Mass-transfer  coefficient  option  control 
variable  (JOPT  =  1  gives  gas-phase  controlled 
coefficient;  JOPT  ■  0  gives  stagnant-drop 
contribution) 

Print  control  variable  (JP  set  equal  to  1 
suppresses  printing  of  individual  raindrop 
data 

Control  variable  for  gas  (=  1,  or  liquid 
(■  2)  based  solubility  calculations 


< 


o 

TABLE  20 

.  (Contd.) 

Symbol 

Units 

M 

r 

O 

N 

P 

atm 

#% 

4> 

PHI 

Q 

gm  moles  /'sec 

R 

cm 

♦  * 

RE 

dimensionless 

REF 

cm 

f  • 

u 

RK 

-1 

sec 

S 

v 

o 

SAVEY 

dimensionless 

SC 

dimensionless 

35 

SIGMA 

cm 

SIGMAY 

cm 

1 

SIGMAZ 

cm 

SIGPHI 

radians 

SIGTHE 

radians 

ft 

SIGY 

cm 

Definition 

Intemal/external  control 

IntPger  for  RUNGE  and  SOLVE 
Number  of  discreet  drop  sizes  in  descretlzed 
spectrum;  also  number  of  simultaneous  equa¬ 
tions  in  RUNGE 

Ambient  pressure 

Internal  computation  variable  in  RUNGE 
Source  strength  of  plume 
Raindrop  radius 

Reynolds  number  for  falling  drop 

Height  above  release  point  where  raindrop 
passes  x  *  0 

First-order  reaction  constant  for  decay  of 
pollutant  in  plume 

(Integer)  status  variable  used  in  subroutine 
RUNGE — counterpart  of  variable  IS  in  sub¬ 
routine  SOLVE 

Mixing  ratio  of  pollutant  in  gas  phase  at 
receptor 

Schiridt  number  for  falling  drop 

Internal  estimate  of  plume  spread  in  z- 
direction — utilized  when  compact  plumes  are 
encountered 

Plume  spread  in  y-direction 
Plume  spread  in  z-direction 
Plume  spread  in  ^-direction 
Plume  spread  in  ©-direction 
Plume  spread  in  y-direction  at  receptor 
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TABLE  20.  (Contd.) 


Symbol 

Units 

Definition 

SIGZ 

cm 

Plume  spread  in  z-direction  at  receptor 

T 

0  K 

Ambient  temperature 

TDZ 

cm 

Compact  plume  grid  spacing 

TEST1 

liter/gm  mole 

Test  variable  for  asymptotic  dilution 
conditions 

TKY 

gm  moles /cm^  sec 

Overall  mass-transfer  coefficient  (cf. 

Eq.  (41)) 

TDZ 

cm 

Grid  spacing  for  compact  plume 

U 

cm/sec 

Mean  wind  velocity 

V 

cm/sec 

Terminal  fall  velocity  of  raindrop  (fall 
in  +z  direction) 

VERT 

cm/sec 

Plume  loft  velocity 

WORATE 

gm  moles /cm 

Downwind  washout  rate 

X 

cm 

Downwind  distance  from  source 

XBUK 

cm 

Downwind  distance  of  receptor  from  source 

XCL 

cm 

X  position  where  drop  falls  across  release 
elevation 

XK 

gm  moles /cm^  sec 

Liquid-phase  mass-transfer  coefficient 

XNT 

cm/sec 

Rainfall  rate 

y 

cm 

Crosswind  distance  from  source 

YAB 

dimensionless 

Mixing  ratio  of  pollutant  in  air  (molts/mole) 

YBUK 

cm 

Crosswind  distance  of  receptor  from  source 

YK 

gm  moles/cm^  sec 

Gas-phase  mass-transfer  coefficient 

Z 

cm 

Distance  above  stack  base 

ZBUK 

cm 

Distance  of  receptor  above  stack  base 
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PROGRAM  EPAEC < INPUT .OUTPUT » TAPE5*INPUT * T APE6*OUTPUT I 
C  PROGRAM  ••EPAEC**  TO  SOLVE  NONLINEAR-NONFEEOBAC*  MODEL  FOR  GAS 
f  SCAVENGING.  PROGRAMMED  IN  FINAL  FORM  MARCH*  1973  BY 
C  J»  M,  WALES  UNOER  CONTRACT  TO  THE  DIVISION  OF  METEOROLOGY. 

C  ENVIRONMENTAL  PROTECTION  AGENCY. 

r 

C  PROGRAM  REAOS  IN  THE  NUMBFN  OF  DISCRETE  DROP  SIZES  IN  THE  RAIN 

C  SPECTRUM  AND  THE  VARIOUS  DECISION  VARIABLES.  IT  THEN  READS  THE 

C  RAIN  DISTRIBUTION  DATA  ANO  (OPTIONALLY)  THE  LIOUIO..PHASE  MASS- 

C  TRANSFER  COEFFICIENTS.  SUBSEQUENTLY  IT  READS  IN  THE  PHYSICAL 

C  PROPERTIES.  THE  BIND  VELOCITY  AND  SOURCE  STRENGTH.  ANO  THE 
C  LOCATION  OF  THE  RECEPTOR.  CALCULATION  THEN  PROCEEDS  AS  DOCUMENTED 
C  THROUGHOUT  THE  PROGRAM,  THIS  CODE  USES  THE  C6S  SYSTEM  OF  UNITS. 

C  EXCEPT  FOR  PRESSURE*  WHICH  IS  IN  ATMOSPHERES.  AND  FOR  SOME 

C  CONCENTRATIONS  EMPLOYED  IN  THE  SOLUBILITY  CALCULATIONS* 
t  WHICH  ARE  GIVEN  IN  MOLES  PER  LITER.  CONVENTIONAL  CONCENTRATION 
C  UNITS  ARE  MOLES  PER  CUBIC  CENTIMETER. 

COMMON  0(20) *F(20) *CGRND(?0) •  SAVEY.N. Jl . J2* J3. J4* JOPT. JP, JENO, 
1DAX,DAY.HEX,XNU»P*T»XNT*SIGTHE*SIGPHI*U.H,O,XBUK,Y0UK,ZBUK. 
20ELTAZ. SIGY.SIGZ. CEO. CAVG.RK.BKG. VERT 
EXTERNAL  V.HPRIME.YAB 

I'M)  READ  (5*320)  N.  Jl  t J2« J3* J4* JOPT  *  JP.  JENO 
IF  (Jl.EO.l)  GO  TO  110 
RF  AO  (5*330)  (D ( I ) *  I *1 *1 0 1 
READ  (5*330)  (F(I). 1*1*10) 

110  IF  (J2.E0.1)  GO  TO  120 

REAO  (5*H0>  DAX.0AY*HEX.XNU«P*T»XNT*RK 
120  IF  (J3.E0.1)  GO  TO  130 

READ  (5.350)  SIGTHE.SlGPHI.UvHfO.VERT.BKG 
130  IF  (J4.EQ.1)  GO  TO  140 

READ  (5.360)  xbuk.ybuk.ibuk.deitay.deltaz 

0YDUM*,0S7266*XBUK 
IF  (OELTAY.LT, 1.)  DELTAY-DYDUM 
140  IF  (J4.EQ.1)  YBUK*0. 

WRITE  (6*230) 

WRITE  (6*240)  DAX.OAY*XNU(HEXtP*T 
WRITE  (6*250)  SIGTHE .SIGPH1 ,U*H*Q 
WRITE  (6*260)  XBUK.YBUKfZHOK 
CCUM*0, 

I  Y*0 

150  call  master 

C  TESTING  FOR  INITIAL  COMPUTATION 

IF  (IY.GT.0)  GO  TO  160 
WRITE  (6*270) 

C  PRINTING  OF  BULK  RESULTS 

160  WRITE  (6.2B0)  SIOY.SIGZfCEG, CAVG 
WRITE  (6*310)  SAVEY 
IF  (JP.EQ.l)  GO  TO  180 
DO  170  1*1 *N 
R*D(I)/2. 

REF*ZBUK-X0UK*V (R) /U-H 
C  PRINTING  OF  INDIVIDUAL  RESULTS 

170  WRITE  (6*290)  D ( I ) «CGRNO ( I ) *REF 
IF ( JEND.EQ.3)  GO  TO  100 

C  TESTING  FOR  WHETHER  RECEPTOR  IS  ON  CENTERL'NE 
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SUBROUTINE  MASTER 

SUBROUTINE  FOR  CALCULATION  OF  *ASMCUT  CONCENTRATIONS  JN  RAIN 
AT  GROUNO  LEVEL  USING  GENERALIZED  NONLINEAR  MODEL.  IT  ACCEPTS 
INPUT  OATA  FROM  MAIN  PROGRAM  THROUGH  THE  COMMON  STATEMENT  ANO 
THEN  PROCEEOS  TO  CALCULATE  RAINDROP  CONCENTRATIONS  Bv  CALLING  Til 
RFOUIREO  SUBROUTINES,  IT  then  CALCULATES  THE  AVERAGE 
CONCENTRATION  IN  THE  RAIN  SAMPLE.  AND  R'-'TU!  NS  CONTROL  TO  THE 
calling  PROGRAM 

COMMON  0(201 ,F(20) , CGRNO ( 2U )  ,  S AVEY ,N. JI . J2. J3, J*. JOPT . JP.dEND . 
IDAXiDAY.HEX,  XNU.P  » T,XNT»Sl(iTHE, SIGPHl  .U.H.O.XBUK  « VBUK  ,  ?RUK, 
2DELTAZ.SlGY,SIGZ.CEQ.CAyG.RK,8KG,VERT 
CALCULATION  OF  STANOARD  DEVIATIONS 
S!GY«X8UK*SIGTHF 
SIGZ«XBl)K. SIGPHl 
DO  100  I > 1  , N 
R«0(I)/2. 

CALCULATION  OF  INOIVIOUAL  CONCENTRATIONS  FOR  RADIUS  »  -’•INOROPS 
CALL  SOLVE  (R.C0BJ.J0PT»DA», DAY. HEX. XNU.P, T.SIGTmp ,SIGPHl ,1  .H.0.X8 
1UK, YflUK.ZBUK, SAVE Y.RK.flKG.vERT. CEO; 

CGRNO ( I ) »COBJ 

CONTINUE 

CAvGmO. 

AAbO  . 

COMPUTATION  OF  AVERAGE  CONLENTRATIO J  BY  DISTRIBUTION  OVER  DROP  ?IZES 
DO  110  I > 1  ,  N 
AAbAA*F ( I ) *0  ( I ) **3 
00  1?0  I  *  1  .  N 

C A vG»CA VG*F ( I)BCGRND(I) *0(1 j *0 

CAvGbCAVG/AA 

RETURN 

FNO 
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FUNCTION  HPRIME (C.HEX.T*P»LGOPT) 

C  SUPPOl  TINE  FOR  CALCULATING  APPARANT  HENRY»S-LAN  CONSTANT  USIN6 
C  JOHNS T ONE -LF PPL A  PARAMETERS. 

HANK»EXP <9, 94-3040, /T 1 /P 
FQCON«EXP(-10.3*1780./TJ 

C  TRANSFER  POINT  FOR  LIQUID-  OR  GAS-PHASE  BASED  CALCULATIONS 
IF (LGOPT.EO.  1 )  60  TO  100 
R»FOCON*HEX 
CO'JH«EQCON*  1  000,-C 

C  TEST  FOR  l OW-C0NCENTRAT10N  CONDITIONS 

TFST1»4.*C0UM/B»*2 

c  ASYMPTOTIC  EXPRESSION  FOR  LOM-CONCENTRAT ION  CONDITIONS 
HPRIMEsiinOO.*  (HANK -HANK •EOCON/BI 
IFiTESTl.LT. .001)  GO  TO  ltu 
C  SOLUTION  OF  TOTAL  SOLUBILITY  EQUATION 

HPRIML*MANK*  < 1 (-R»SORT («••?. 4, aCDuM) )/ (2000 . *C )>• 1 000 , 

GO  TO  110 

C  TEST  TOR  LOW-CONCENTRATION  CONDITIONS 

100  DUMMVsHEX«HEX*HANK/  (4*FQC0N«C) 

IF(OUMMY.LT.1000.)  GO  TO  105 

C  ASYMPTOTIC  EXPRESSION  FOR  LOw-CONCENTRaT ION  CONDITIONS 
HPR IME» 1 000«*HANK»hEX/( NE  *  *EQCON ) 

GO  TO  110 

C  SOLUTION  OF  TOTAL  SOLUBILITY  EOUATION 

105  CTEST* (C/HANK* (-HEX* SORT (HtX**2»4*C*EQCON/HANK) )/2)/1000. 
HPRIME-C/CTFST 

C  RETURN  VALUE  OF  HPRIME  HAS  UNITS  OF  CENTIMETERS  CUBED  PER  MOLE 
110  RETURN 
ENn 


|  The  above  subroutine  is  that  employed  for  the  prediction  of  SO2  washout 

from  power  plant  plumes.  The  sensitivity  of  SO2  solubility  to  pH  required 
the  incorporation  of  pH  data  and  the  SO2-H2O  equilibrium  constant  to  the 
calculation  of  Henry's  Law  constant,  HPRIME.  Application  of  the  EPAEC 
P  model,  to  the  scavenging  of  ethylacetoacetate  and  diethylamine  involved 

input  of  H'  values  calculated  in  Appendix  A. 
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SUBROUTINE  SOLVE  <R,COBJ.JOPr,DAX.PAV.Hex>XNU.P.T.SlGTME.SlGPMt.U. 
i H  tQf  XBUK »  Y.  ZBUK  * SA VE Y ,RK . tt*G . VERT  »C£Q ) 

C  SUBROUTINE  FOR  CALCULATION  OF  INDIVIDUAL  RAINDROP  CONCENTRATIONS* 

C  PROGRAM  SETS  UP  APPROPRIATt  COMPUTATION  GRID  FOR  NUMERICAL 
C  SOLUTION  OF  THE  FIRST«ORDfN*  ORDINARY  DIFFERENTIAL  EQUATION 
C  DESCRIBING  DROP  RESPONSE.  IT  THEN  CALLS  THE  RUNGE-KUTTa 
C  ALGORITHM  REPEATEDLY.  PROGRESSING  FROM  THE  TOP  OF  THE 
C  COMPUTATION  GRIO  TO  GROUND  LEVEL. 

DIMENSION  C  1 1 >  «F  1 1 ) 

c  initialization  of  variables 

ICP-0 

M»0 

IGRNDaO 

C(l)-0 

I«1 

C  BYPASS  TEST  FOR  EQUILIBRIUM  SCAVENGING  CONDITIONS 

HSTAR«H*V£RT«XBUK/U 

sigmay«sigthe*xbuk 

SIGMAZ«SIGPHI»XBUK 

SAVEY>YAB(T.P.U»HSTAR»Q«XBUK.Y.ZBUK»SIGTHE.SIGPHItRK*BKG) 

HTEST-HPR IME ( SAVE Y, HEX, T.P.l) 

IF (HTEST.LE.O.)  GO  TO  IOO 
CTEST-SAVEY/HTEST 

CALL  MTC(R.OAX.DAY*XNU.T.P«dOPT,CTEST,HTEST.TKY» 
GR«1S.9cTKY«HTEST*U*SI6MAV*SIGMAZ**2*EXP ( »5*.5*Y*Y/SIGHAY**2) 
6Ra«GR/(V(R|*R> 

GPOUP»GR*SAVEY/  (8i?,*0*T/P) 

CORJ-CTEST 

CEO.CTEST 

IF (GROUP. GT . IS. )  GO  TO  190 
C  CALCULATION  OF  INITIAL  RAINOROP  CONCENTRATION 
C II » •BKG/HPRIME lBKQ.HEXtT.Pf 1 ) 

C  TEST  FOR  PLUME  UNDERCUT 

100  REF*-XBUK«V<R)/U«ZRUK-H 
Z*PEF*H 

IF  (REF. GT. 0.1  GO  TO  120 
C  SPACING  FOR  UNOERCljT  GRID 

OE  ?AZ»»XBUK*V(R)/(U*30JJ 
60  TO  130 

C  SETTING  OF  NORMAL  GRIO  SPACING 

120  XCL«XBUK«U*(H-ZBUK)/V(R| 

SIGMA«SIGPHI*XCL 
HSTAR«H*VERT*XCL/U 
Z1«HSTAR«3.*SIGMA 
IF  (Zl.LT.Z)  Z«ZJ 
OELTAZaZ/60. 

C  TEST  FOR  COMPACT  PLUME 

0ZTST«S1GMA/A. 

IF (SIBPHI.GT..S)  OZTST«DZT ST/2, 

IF (DZTST.GT .DELTAZ)  GO  TO  130 
C  SETTING  OF  COMPACT  GRIO  SPACING 

TOZ*OZTST 
ICPal 

C  START  OF  NUMERICAL  INTEGRATION  LOOP 

130  IF  (ICP.NE.l .0R.I.GT.2S)  GO  TO  140 
OZ«-TOZ 


r 
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z 


z 


t: 


r. 


60  TO  ISO 
140  DZ»«DELTAZ 

C  INITIATION  OF  RUN6E-KUTTA  ALGORITHM 

iso  call  bunoe  ii*c»f»z*oz«is*mi 

IF<C<i>.LT.O.)  Cd)«0, 
ir  (IS.NE.l)  60  TO  1*0 
X«X8UK«U« (Z-Z6UK ) /V (R) 

HSTAR«H*VERT«X/U 

C  CALCULATION  OF  MASS-TRANSFER  COEFFICIENT 

CALL  MTC  <B«0AX.0Av.XNU«T.F.J0BT#C(1)  .HEX»TKV» 

C  DETERMINATION  of  FIBST  DERIVATIVE  FOB  B-K  ALGORITHM 

F(1)»(3.*TKY/|R*V  CR)  )>•  »Y*mT.P.UtHSTAR.0»X#y.Z#SI6TMEtSI6PHl.RKiH 
1K0|-HRB1ME(C<1) *MEX*TiF«2t *C(ln 
60  TO  150 

C  TEST  FOB  APPROACH  TO  GROUNO 

160  ZTEST»Z*OZ-ZRUK 

tr<ZTEST.LT.100.J  60  TO  170 
I  •  I  ♦  l 
60  TO  130 

170  IF  (IGRND.Ea.il  60  TO  100 

I6RN0*1 
OZ-ZBUK-Z 
60  TO  ISO 
160  C0BJ»C(1 ) 

190  RETURN 

ENO 


SUBROUTINE  MTC  (Rf DAX«OAYt XNU*T»Pf JOPT .CSET • HEX* TRY! 

C  SUBROUTINE  FOR  MASS-TRANSFfcR  COEFFICIENT  CALCULATION. 

C  OAS  COEFFICIENT  BA$EO  ON  FM0ESSLIN6  EQUATION.  LIQUID 
C  COEFFICIENT  IS  BASEO  ON  CONTINUITY  EOUATJON  SOLUTION  FOR 
C  RESPONSE  TO  RAMP  CONCENTRATION  FORCING  FUNCTION. 

C 

C  CALCULATION  OF  3AS.PHASE  COEFFICIENT 

RE»«t.*R*V (R) /XNU 
SC«XNU/DAV 

YK* ( l.».3*RE**.5*SC«*«333) •DAY*P/(R*T«6?.057) 

TKY-VK 

IF  (JOPT.EO.l)  60  TO  100 

C  CALCULATION  OF  LIQUlO-PMASt  COEFFICIENT 

XKa.Z77B«0AX/R 

TRV«l/(HPRIME|CSET»HEX»T«PtZ)/XK»l/YK) 

100  RETURN 

ENO 
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SUBROUTINE  RUNOE  <N,YtF*X*H*$*M) 

C  SUBROUTINE  TO  SOLVE  DIFFERENTIAL  EQUATION  DESCRIBING  CONCENTRATION 

C  RESPONSE  OF  A  FALLING  ORO*.  RUNGE-KUTTA  ALGORITHM  1$  ADAPTED  FROM 

C  GENERAL  VERSION  GIVEN  BY  CARNAHAN*  LUTHER*  ANO  WILKES* 

DIMENSION  SAVEV ( 1 ) *PHI  (I) *F (1) *Y(1) 

INTEGER  S 
mbm*i 

GO  TO  (100*110*130*150*170)*  M 
100  Sal 

GO  TO  190 

110  00  120  J«1*N 

SAVEY(J)«V(J) 

PHI ( J)aF(J) 

120  Y(J)aSAVEY<J)*.5«H«F(J) 

XaX«.S«H 

S«1 

GO  TO  190 

130  00  1*0  Jal.N 

PHI(J)apHl(J)«2*aF|J) 

1*0  V ( a$AVEV ( J) **S*H«F ( J) 

Sal 

GO  TO  190 

150  00  160  Jal,N 

PHI(J)apHI(J)«2,aF(J> 

160  V (J)aSAVEY ( J) *H«F ( J) 

XaX«,5aH 

Sal 

GO  TO  190 

170  DO  160  Jal,N 

160  Y (J)aSAVEV ( J) ♦ (PHI (J) *F |JI ) aH/6* 

NaO 

Sa? 

190  CONTINUE 
RETURN 
END 


